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REFERENCE PUBLICATION 


BASIC MECHANICS OF LAMINATED COMPOSITE PLATES 


I. INTRODUCTION 


A. Intent and Scope 

This report is intended only to be used as a quick reference guide on the mechanics of 
continuous fiber-reinforced laminates. By continuous fiber-reinforced laminates, the following is 
assumed: 

(1) The material to be examined is made up of one or more plies (layers), each ply consisting 
of fibers that are all uniformly parallel and continuous across the material. The plies do not have to be 
of the same thickness or the same material. 

(2) The material to be examined is in a state of plane stress, i.e., the stresses and strains in 
the through-the-thickness direction are ignored. 

(3) The thickness dimension is much smaller than the length and width dimensions. 

An attempt is made in this report to develop a practical guide that can be easily referenced by 
the engineer who is not familiar with composite materials, or to aid those who have seen this subject 
matter before. 

The scope of the report will be limited to the elastic response of the above-mentioned class of 
material. Strength-of-laminated composites will not be covered. General composite material 
mechanics and strength are developed in more detail in texts such as Jones 1 and Halpin. 2 It is 
assumed that the reader has a general knowledge of elastic stress-strain behavior. 


B. Terminology and Notation 

Some terminology important to composite materials follows: 

Isotropic— Possessing the same mechanical properties in all directions. Composite laminates 
are never isotropic. 

Laminate — A material consisting of layers (laminae) bonded together. 

Transversely Isotropic — Possessing one plane that has the same mechanical properties at 
any direction in that plane, i.e., the laminate will have the same stress-strain behavior at any 
direction in the plane of the material (sometimes called quasi-isotropic). 

Orthotropic — A material that has different mechanical properties in three mutually 
perpendicular planes. Note that the properties of the material are direction specific in this case. All 
unidirectional laminae are individually orthotropic. Most laminated composites fall into this category. 



Homogeneous — Material properties do not change from point to point within the material. 
Since filamentary composites consist of at least two distinct phases (fiber and matrix), laminated 
composites are never truly homogeneous, although on a macroscopic scale when discussing the 
linear elastic response (no damage) of laminated composites the material may be generalized as 
homogeneous. This assumption is termed “smearing” of fiber and matrix. In reality all composites 
are heterogeneous. 

Principal Material Directions — Directions parallel and perpendicular to the fibers in a lamina. 
Note that these directions are not necessarily the directions of principal stress as defined by 

continuum mechanics. 

Balanced Laminate — For each +0 ply in the laminate there is an equally thick -6 ply in the 
laminate. This does not apply to 0° and 90° plies. 

Symmetric Laminate — The plies of the laminate are a mirror image about the geometrical 

midplane. 

Angle Ply Laminate — Containing plies oriented at angle(s) other than 0° or 90°. 

The notation used throughout this report denotes the directions parallel, perpendicular, and 
through-the-thickness to the aligned fibers of a ply of material as 1,2, and 3, respectively. The 
coordinates of the laminate are denoted by x,y (in-plane) and through-the-thickness z- Some texts 
and laminate analysis computer software have reverse notation, which is exactly opposite to that 
presented here. Some users prefer to use a primed and unprimed notation. Thus, it is important to 
note which notation is being used. 

Since a state-of-plane stress is assumed, the term “plate” will be interchangeable with 
material and composite laminate. Please note that in actual practice, the laminate does not need to 
be a plate, but may be a shell or other shape such that the material is considered “thin.” 


C. Summary of Sections 

In section n, the behavior of an individual ply or lamina is considered. This behavior is the 
“building block” upon which laminated plate theory is based. The important concept of coupling, 
unique to anisotropic materials, will be introduced in this section. 

The main body of this report is considered in section III. It develops the relationship between 
loads and deflections for a composite laminate. Most information on the elastic response of 
laminated composites can be calculated from the equations given in this section. The constitutive 
equations that govern the laminate load/deflection behavior are given in this section. 

For those who have not dealt with laminates, section IV shows how stacking sequences of 
the plies that make up the laminate are denoted. This section simply presents a standard “code” by 
which a laminate can be described. Many will already be familiar with this material. 

Section V contains derivations of the most sought after numbers when dealing with the 
elastic response of these materials, the engineering constants. Many texts on laminated composites 
omit this part, or only present it for special types of laminates that make the calculation relatively 
simple. The engineering constants for any laminate can be found from the equations in this section. 


2 



Section VI introduces the effects of temperature and moisture on the strain of a composite 
laminate. These effects are often neglected, but are very important in determining the stresses and 
strains within each ply of the laminate. Determining these ply stresses and strains for symmetric 
laminates is presented in section VII and in section VIII for unsymmetric laminates. These sections 
provide information necessary to study the strength of composite laminates. 

Section IX is a summary of the most important equations presented in this paper. 


II. GENERALIZED HOOKE’S LAW FOR NONISOTROPIC MATERIALS 


A. Normal Stress and Strain, Uniaxially Applied Force 

Normal stress is defined as the force per unit area acting perpendicular to the surface of the 
area. The corresponding strain is defined as the elongation (or stretch) per unit length of material in 
the direction of the applied force. For isotropic materials, the relationship between stress and strain 
is independent of the direction of force, thus only one elastic constant (Young’s modulus) is required 
to describe the stress-strain relationship for a uniaxially applied force. For a nonisotropic material, 
at least two elastic constants are needed to describe the stress-strain behavior of the material. 

Figure 1 is a schematic of an isotropic and a unidirectional fiber-reinforced material. The 
stiffness of the isotropic plate can be described by one value, the modulus, E, of the material, 
regardless of direction of load. The stiffness of the orthotropic plate must be described by two values, 
one along the longitudinal direction of the fibers, commonly referred to as El, and one transverse to 
the direction of fibers, usually denoted by Ej- Subscripts 1 and 2 will be used such that El- E \ and 
E t = E 2 . Thus, indices must be added to the stress, strain, and modulus values to describe the 
direction of the applied force. For example, for an isotropic material, the stress/strain relationship is 
written: 


o=E e 


( 1 ) 


Isotropic Plate 


Orthotropic Plate 


1 


k 



stiffness in 1 -direction = stiffness in 2-direction 
= stiffness in any direction 


Reinforcing fibers 



stiffness in 1 -direction » stiffness in 
2-direction * stiffness in other directions 


Figure 1. Difference between an isotropic and an orthotropic plate. 



For the orthotropic system, the direction must be specified. For example: 

<7j = £ ! £j or (T 2 = £ 2 e 2 • (2) 

If the applied load acts either parallel or perpendicular to the fibers, then the plate is considered 

specially orthotropic. 


B. Stress and Strain, Plane Stress for Specially Orthotropic Plates 

The previous section dealt with an extremely simple type of stress state, uniaxial. In general, 
plates will experience stresses in more than one direction within the plane. This is referred to as 
plane stress. In addition, Poisson’s ratio now becomes important. Poisson’s ratio is the ratio of the 
strain perpendicular to a given loading direction, to the strain parallel to this given loading direction: 


Poisson’s ratio = v 


12 ' 


c 2 

— or 


°1 

For loading along 
the fibers 



For loading 

perpendicular to the fibers 


(3) 


The strain components are now stretch due to an applied force, minus the contraction of 
Poisson’s effect due to another force perpendicular to this applied force. Thus: 


Using equation (2): 



V 21 f 2 


and e 2 = F _V 12 £ l • 



and e 2 




(4a) 


(4b) 


Shear forces can also be present. Shear stress and shear strain are related by a constant, like 
the normal stresses and strains. This constant is called the shear modulus and is usually denoted by 

G. Thus: 


T 12=Tl2 G 12 • (5) 

Where tn is the shear stress (the 1 and the 2 indices indicating shear in the 1-2 plane), and J \2 is 
the shear strain. Figure 2 gives a definition of shear strain. 

Since it is known that a relationship exists between Poisson’s ratios and the modulii in each 
of the two axes directions, namely: 


v 21 E l = v n E 2 , 


(6) 




(b) 

Figure 2. Definition of shearing strains. 
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Equations (4b) and (5) can be written in matrix form as: 


e l 


611612 0 


^1 

£ 2 

= 

s 12 22 0 


<*2 

/l2 


0 0 6 66 


7 12 


where, 

S„ = l/£i S 22 =\IE 2 

S l2 = -v l7 /E l = -v 2] /E 2 5 66 =1/G 12 


(7) 


( 8 ) 


Note that at the 3,3 position in this 3x3 matrix (called the compliance matrix), the subscripts 
are 6,6. This evolves from a detailed treatment of arriving at a constitutive equation for an orthotropic 
material from an anisotropic one. 

By inverting the compliance matrix, one can get stress as a function of strain. This turns out 

to be: 


*1 


011012 0 


£ 1 

°2 

= 

012 022 0 


£ 2 

T 12 


0 0 q 66 


n 2 


where: 


0 


it 


012“ 


1 " V 12 V 21 

1 “ V 12 V 21 1 - V 12 V 21 


022 “ 


l ~ V l2 V 21 


066-^12 


The Q’s are referred to as the reduced stiffnesses and the matrix is abbreviated as [£>]. 


(9) 


( 10 ) 


C. Stress and Strain, Plane Stress for Generally Orthotropic Plates 

Now suppose that the unidirectional lamina in figure 1 is loaded at some angle other than 0° 
or 90°. The lamina is now referred to as generally orthotropic, since, in general, the loading direction 
does not coincide with the principal material directions. The stresses and strains must now be 
transformed into coordinates that do coincide with the principal material directions. This can be 
accomplished using the free-body diagram in figure 3. From free body diagram (a), and summing 
forces in the 1-direction: 
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2 


1 


Gx 


T12 y 



Area of face = dA 
<71 



Area of face = dA 



ay 

(b) 

Figure 3. Generally orthotropic lamina. 

= 0 = o x dA-aJdA cos tfjcos 6 - cr y [dA sin sin 0 
-<J. dA cos dj sin sin cos 6 

From the free body diagram of figure 3(b) and summing forces in the 2-direction: 

£F 2 = 0 = <j 2 dA-cJdA sin sin 6 -cr y {dA cos fljcos 6 
+t 'xy[dA cos fljsin 6 + sin #jcos 6 
From free body diagram of figure 3(b) and summing forces in the 1-direction: 

2 F 2 = 0 = T n dA+crJdA sin cos 6 -crJdA cos sin 6 
-r J^dA cos cos 6 +z x jdA sin 6?jsin 6 


►x 


( 11 ) 


(12) 


Simplifying equations (11), (12), and (13); 


o, = o cos 2 9+o sin 2 Q+ 2r sin 9 cos 9 , 

* x y *y 

O-, — o sin 2 9+o cos 2 9- 2t sin 0cos 0 , 

2 * y .xy 

Tj 2 =- ct. sin # cos tf+cr^sin 9 cos 9 +t^cos 2 9- sin 2 e) . 
Equation (14) can be written in matrix form as; 


1 


2 ^ 

. 2 „ 






cos 9 

sin 9 

2 sin 9 cos 9 


°x 

°2 

= 

■ 2 o 

sin 9 

2 

cos 9 

-2 sin 9 cos 9 


°y 

T 12 


- sin 9 cos 9 

sin 9 cos 9 

( cos 2 9- sin 2 £?j 




( 14 ) 


(15) 


The 3x3 matrix in equation (15) is called the transformation matrix and is denoted by [T). The same 
matrix is used to transform strains. Note that the tensorial shear strain must be used, not the 
engineering shear strain, when transforming strains. This arises from the geometrical considerations 
that the amount of shear must be equivalent with respect to both the x- and y-axes, since these 
axes will be transformed into new ones (fig. 2(b)). 


If it is desired to transform from the 1-2 coordinate system to the x-y coordinate system, the 
inverse of [ 7] must be found. It is given by: 


Thus; 



2 ^ 
cos 9 

sin 9 

sin 9 cos 9 


■ 2 

sin 9 
2 

cos 9 
- sin 9 cos 9 


-2 sin 9 cos 9 
2 sin 9 cos 9 
| cos 2 9- sin 2 9j 


°\ 


a 


<7 


a i 


X 


X 


CT 2 

= m 

°y 

and 

°y 

= m"' 

°2 

T 12 


> 


> 


T 12 


(16) 


( 17 ) 


Similarly for strain; 


e l 






e i 

£ 2 

= m 


and 


=m" 1 

£ 2 

£ 12 




£ xy 


£ \2 


Putting equation (9) into the second part of equation in (17): 



[r]-‘[e] 


£ 1 


100 

e i 

£ 2 

=[7-]- I [e] 

0 1 0 

£ 2 

y\ 2 

002 

£ 12 


Now putting the first equation of equation (18) into equation (19): 


( 18 ) 


°x 


100 



°y 


0 1 0 
002 

M 

£ xy 


Defining a new matrix called the lamina stiffness matrix (sometimes called “(2-Bar”) as: 


( 20 ) 


and letting: 


[sHrHe] 


1 0 0 
0 1 0 
0 0 2 


[T] ■ 


( 21 ) 


m = cos 6 , 


the components are: 


n = sin 6 , 


2 it = 2iim 4 +2((2 1 2+2G 66 )m 2 n 2 +2 22 n 4 , 

2 12 = (Q h+022- 4 266) mn+Q 12 (m 4 +« 4 ) , 

222 = <2n« 4 +2(2i2+2G 66 ) rnn 2 +Q 22 rn , 

2 16 = (2ll~2l2~2266) m n + {Ql2-Q22 + 2Q66) m,1 ’ 


( 22 ) 



2 26 = ( 0 11 -Q 12 _ 2 266 ) n3/n + [Q 12 ~ 022 + 2066 ) nm > 
066 = [Qu + Q22~2Q 12 — ^S 6e) m n + 06e( m +n ) • 


(22) (cont.) 


Note that if 0is any angle other than zero, there will be nonzero 0 16 and 0 26 terms - Putting this 

into equation (20): 



011 

012 
016 


012 2016 
022 2026 
026 2066 



a x 


011 012 016 

<y 

= 

012 0 22 026 

.v 


016 026 066 



(23) 


it can be seen that a shear strain will produce normal stresses, and normal strains will contribute to 
a shear stress. This is referred to as extension-shear coupling and will take place in a lamina that is 
loaded at an angle to the fibers (other than 0° and 90°). That is, there will be coupling if the 0 16 
and/or 0 2 6 terms in the lamina stiffness matrix are nonzero. 


D. Invariant Stiffnesses 
The “0-Bar” terms can be written as: 


0 11 = [/ 1 +[/ 2 cos(20)+[/ 3 cos M. 
022 = ^l~^2 cos (2e)+u 3 COS M. 

012= U4-U3 COS (40) , 

0 66 = £/ 5 -f/ 3 cos(40) , 

G 16 = |t/ 2 sin(20)+t/ 3 sin(4e) , 
0 26 = |f/ 2 sin(20)-t/ 3 sin(4^) , 
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where. 


3 11 

^l-^(Qn + Q22) + ^Ql2 + 2^66 » 

^2 = ^ (2ll _ 222) ’ 

= — (2ll+222)-^ 2 12~~ 266 > 

^4 = ^ (2ll+222) + ^ 2 12- 2 266 » 

1 11 

^5 = ^(211+222)-^ 2l2 + 2 266 • 


(24b) 


Note that only C /2 and C /3 are coefficients to the sine or cosine terms in equation (24a). This implies 
that when calculating the <2-Bar values, U 1 , C/ 4 , and C /5 are independent or invariant to the ply 
orientation 6. This concept of “invariant” quantities can make some calculations easier. This paper 
will not go into detail on this subject since only the basics are being presented. 


HI. MECHANICS OF LAMINATED COMPOSITES 


A. Assumptions 

The following assumptions are made for the remainder of this paper: 

(1) The laminate thickness is very small compared to its other dimensions. 

(2) The lamina (layers) of the laminate are perfectly bonded. 

(3) Lines perpendicular to the surface of the laminate remain straight and perpendicular to 
the surface after deformation. 

(4) The laminae and laminate are linear elastic. 

(5) The through-the-thickness stresses and strains are negligible. 

These assumptions are good ones as long as the laminate is not damaged and undergoes 
small deflections. 


B. Definitions of Strains and Displacements 

A displacement of the plate in the ^-direction is designated as u. For the y-direction, it is 
designated as v and for the z-direction w. Figure 4 shows these displacements. The strains are now 
defined as: 
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y 



|s Deformed 

shape 



1 

1 

Original 

shape 

1 

1 

1 

1 

1 

1 

1 

1 


u U- 




Normal displacement 

z 





Figure 4. Displacements of a plate. 


du 

dx 







(25a) 


The slope of the plate if it is bending is given as: 


— along the x-direction , 
ax 


— along the ^-direction . 

dy 


(25b) 


The total in-plane displacement at any point in the plate is the sum of the normal displacements plus 
the displacements introduced by bending. Denoting the displacements of the midplane of the plate for 
the x and y directions as u 0 and v 0 respectively, with the help of figure 5 the total displacements are: 


u = u 0 -z 


dw 

dx 


v = v 0 -z 


dw 

dy 


(26) 


Note that for figure 5, x can be replaced by y when u is replaced by v (i.e., the view could be from any 
side of the plate). It is assumed that there is no strain in the thickness direction, only a displace- 



z 


Undeformed edge of plate 


Top surface of plate 


: : ( 

+z 

— 

-z 

i — — 



x 


0 = Slope = 


9w 

9x 


For small 0 


Displacement along z-direction 

due to bending is sin9 times z. 

Since 0 is small, sin0 = 0. 

Therefore the displacement is 
- z 0 where negative is compression and 

positive is tension 


Figure 5. Total displacements in a plate. 


From equations (25) and (26): 


Defining: 


du 

Bu 0 

-,2 
d w 

£ = = 


-Z 0 

* dx 

dx 

dx 

dv 

Bv a 

d 2 w 


£y By By Z dy 2 


->2 
d w 


du dv Bu 0 dv 0 
Yxy = ~dy + d~x~~dy + Z &dy 


Bu O O Bv 0 0 

as e ; — as e ; and 

Bx 0 By , 0 y 


Bu 0 dv 0 0 

+ as r 

By a dx 0 r v 


to be the midplane strains and defining: 

as K x 


->2 
d w 


dx 


n 2 
d w 

— - as K ; and -2 

By 2 


■,2 
d w 

as K„, 

dxdy 39 


( 27 ) 


to be the plate curvatures will make notation easier. 
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on (27) can now be written in matrix form as: 




O 



e x 


£ x 


K x 

e 


o 

£ 

+ Z 

K v 

y 


y 

y 

r iy. 


o 

7 

•xy 


Kxy 


(30 


figure 6, it can be seen that the plate curvature K x or K y is the rate of change of slope of 
•late in either the x- or ^-direction, respectively. The plate curvature term Kjy is the 
nding in the x-direction along the y-axis (i.e., twisting). 


k 

k 




From equation ( 23 ), the stresses in each ply of the laminate can be determined with equation 


( 30 ): 


*y 


Qu <2 12 (2i6 
Qn 0.22 2 26 
016 2 26 <266 



+ z 


0 it Q12O16 
0 12 0 22 Q 26 
016 026 066 


K x 


K 


xy 


(30a) 


C. Definitions of Stress and Moment Resultants 

Since the stress in each ply varies through the thickness of the laminate, it will be convenient 
to define stresses in terms of equivalent forces acting at the middle surface. Referring to figure 7, it 
can be seen that the stresses acting on an edge can be broken into increments and summed. The 
resulting integral is defined as the stress resultant and is denoted by N;, where the i subscript 
denotes direction. This stress resultant has units of force per length and acts in the same direction 



Total force in x-direction = 2 Gx(dz)(y) 



Figure 7 . Definition of stress resultant. 
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as the stress state it represents. Figure 7 could also be drawn for the y-direction stress and shear 
stress. The three stress resultants are therefore: 


rm 


N = 

c^dz , 

x J-m 


rm 



adz , 

y J-m 

y 

rm 


Af = 

t dz ■ 

39 J-m 

xy 


(31) 


As can be seen from figure 7, the stress acting on an edge produces a moment about the midplane. 
The force is a x ( dz)(y ) as denoted in figure 7. The moment arm is at a distance z from the midplane. 
Following the same procedure as for the stress resultants, the moment resultants can be defined as: 

rh/2 

M x = ojzdz , 

J-m x 

rhJl 

M = azdz , (32) 

y J-m y 

rh/2 

M= T zdz . 

39 J-m ** 


These moment resultants have units of torque per unit length. 

The directions for all of the stress and moment resultants are shown in figure 8 for clarity. The 
double-headed arrow indicates torque in a direction determined by the right-hand-rule (i.e., point the 
thumb of your right hand in the direction of the double-headed arrows and the direction of rotation of 
the torque is in the direction that your four fingers are pointing). Note that M x and M y will cause the 
plate to bend and M xy will cause the plate to twist. 

As an example of the relationship between stress and stress resultants, if a tensile test 
specimen is 2.54-cm (1-in) wide and 2-mm (0.08-in) thick, and is pulled on with a force of 4,500 N 
(1,000 lb), then the average stress on the cross section is: 




4,500 N 


= 88.6 MPa = 12,500 


lb 


(o.0254 m) (0.002 m) 

The stress resultant would be the average stress multiplied by the specimen thickness: 

N x = (88.6 MPa) (o.002 m) = 177.2 ^ = 1 ,000 ^ 


(33) 


(34) 
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Stress Resultants 



Figure 8. Direction of stress and moment resultants. 

Since 1,000 lb was applied over 1 inch of specimen length, the expected result of 1,000 Ib/in is 
obtained. Each individual ply of the specimen may have a stress other than 88.6 MPa (12,500 
but the average stress will be 88.6 MPa (12,500 lb/in 2 ). 


D. Constitutive Equations for a Laminate 


Putting equation (31) in matrix form: 



and putting equation (32) in matrix form also: 


M x 


rhJ 2 

a x 

M v 



a 

y 


/-A/2 

y 

, s 





zdz 


(36) 


The integrals in equations (35) and (36) must be performed over each ply and then summed, 
since discontinuities in stresses can occur at ply interfaces. Using the schematic of a laminate in 
figure 9, equations (35) and (36) must be written as: 


and. 



dz 


zdz 


(37) 


(38) 


Now equation (30) can be substituted into equation (23), which can then be substituted into 

equations (37) and (38) to give: 




t 


r 

1 



- 


0 

N x 



l2n Q 12 Gw 



Ny 

- k h\ 

rh 

r h k 
^ h k- 1 

Q 12 G22 Q 26 


0 

> 

-JT 


Q 16 Q 26 2 66 

k 

0 



V 


L 


+ 



1 



\ 

fill fi 12 file 


K x 


fil2 fi22 fi26 


Ky 

zdz \ 

2 16 226 066 

k 

Kxy 



( 39 ) 


and, 



Ply k+1 

Geometric 
roidphne- 


Ply k 


Ply k-1 


Ply 2 


Ply 1 


Note that ply k and ply k+1 are the same lamana (layer), but are separated into 
two plies by the geometric midplane. 

Figure 9. Cross section of a laminate. 


M x 


M. 


xy 


= x f 

k = 1 


fill fil2 fil6 
fi 12 Q22 Q26 
fi 16 Q26 Q66 


J k 


y 

0 

Yv 


zdz 


+ 


r h ic 


fill fi 12 fi 16 


K x 

2 1 

fi 12 fi22 026 


Ky 

z dz\ 

fi 16 fi 26 fi66 

k 

K*y 

1 


( 40 ) 
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Since the middle surface strains and curvatures (the £°’s and K’s) are not a function of z (because 
these values are always at the middle surface z = 0), they need not be included in the integration. 
Also, the laminate stiffness matrix is constant for a given ply, so it too will be a constant over the 
integration of a lamina thickness. Pulling these constants to the front of the integral in equations (39) 

and (40) gives: 


and 


N x 

K 




xy 


n 

-x 

k= 1 


2 11 Qn Ql6 
Q12Q22Q26 
2 16 2 26 066 



r _ _ 


r -I 


2ll 2l2 2l6 


K x 

f h k , 1 

2l2 222 226 


Ky 

zdz ) 
Vi / 

2l6 2 26 266 

k 

Kxy 


M x 

M y 

Mxy 



| Gllfil2Gl6 
\ 2 12 Q-22 2 26 
1 2 16 226 266 



r _ _ _ 1 


r -I 

\ 

2 11 2 12 2 16 


K x 

C h k 2,1 

Q 12 Q 22 Q 26 


K y 

z dz ) 

Vi | 

Ql6 226 266 

k 

Kxy 


(41) 


(42) 


Performing the simple integrations gives: 
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xy 



2n 2i 2 2i6 
2 12 Q 22 Q 26 
2 16 226 266 


Jjk 


y 

0 


[ h k~ h k- 1) 



r— — — 1 





2n 2 12 2 16 


K x 

1(2 2)1 

+ 

2 12 2 22 2 26 


K y 



216 226 266 

k 

Kxy 

J 


(43) 


and 
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My 
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( 44 ) 


Since the middle surface strains and curvatures are not a part of the summations, the laminate 
stiffness matrix and the h k terms can be combined to form new matrices. From equations (43) and 
(44), these can be defined as: 


n 

A ij~ ^ ^£-l) » 



In matrix form, the constitutive equations can easily be written as: 


(45) 


(46) 


(47) 


N x 


An ^ 12^16 

B n #12 

N y 


A 12 A 22 A 26 

B 12 B 22 



A 16 A 26 A 66 

5 16 B 26 

M x 


B 11 B 12 B 16 

^11 ^12 

M y 


B 12 B 22 B 26 

D 12 D 22 



B 16 B 26 B 66 

D 16 D 26 


B 16 
B 26 
B 66 


£ 

e 

£ 


xy 


D 


16 


^26 

D 66 


K 


xy 


(48) 


Written in contracted form, equation (48) becomes: 
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A 

1 s’ 


O 1 
CO, 

.M _ 


B 

1 D 


[K\ 
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This can be partially inverted to give: 




A* 

B*' 


' N ' 

M 


c* 

D* 


. K _ 


where. 


[a'UaF, 

[«*]— MT*] • 

[ci-wur 1 . 

M=[b]-[s][a]' 1 [s] . 


(50) 


(51) 


The fully inverted form is given by: 


where. 


A 

K x 

K y 


A 11 A 12 A 16 

/ / / 

A 12 A 22 A 26 
A 16 A 26 A 66 


C 11 ^12 ^16 
/ / / 

C 12 C 22 C 26 
C 16 C 26 ^66 


B 11 B 12 B 16 

B\2 B 22 B 26 
B 16 5 26 5 66 


Du D l 2 D 16 

t f / 

D n E>22 D 26 

D 16 D 26 D 66 


N x 

N y 

Kxy 

M x 

My 

M 


xy 




M-wr . 


[c\-[dT[c] , 


(52) 


(52a) 


[dUd] 


-1 


The fully inverted form is the most often used form of the laminate constitutive equations. 



For symmetric laminates (laminates that are configured such that the geometric midplane is a 
mirror image of the ply configurations above and below the midplane), the geometric midplane is also 
the neutral plane of the plate, and the [5] matrix will have all elements equal to zero (as will be 
shown later). However, if the laminate is unsymmetric, i.e., if the plies near the bottom of the plate 
are much stiffer in the x-direction, then the geometric midplane will not be the neutral plane of the 
plate; and the neutral plane will be closer to the bottom of the plate for x-direction bending as shown 
in figure 10. This is accounted for in the constitutive equations, since the [B] matrix will have some 
nonzero elements (as will be shown later), implying that a bending strain (plate curvature) will 
cause a midplane strain as depicted in figure 10. Likewise, a midplane strain will cause a bending 
moment. A method to find the neutral axis of the plate will be discussed in a later section about 
stresses within the plies of a laminate. 


Undeformed edge of plate 

Top surface of plate 



+z 


-z 

Very stiff plies 




Figure 10. Displacements in an unsymmetrical plate. 

E. Physical Meanings of the [A], [B], and [D] Matrices 
Recalling the definitions of the [A], [B], and [D] matrices. 


A ij=lL[Qij\ k [ h k-hk-i) > 

/fc= i 


n 


l 



( 45 ) 


( 46 ) 


( 47 ) 
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and referring to figure 9, it can be seen that the last term in equation (45) is the kth lamina thickness 
which will be denoted by t k . Thus, equation (45) can be written as: 



(53) 


This matrix is called the extensional stiffness matrix. From the constitutive equation (48), it 
can be seen that these terms relate the normal stresses and strains (much like the modulii of 
elasticity), except for the A 16 and A 2 6 terms which relate shear strains to normal stresses and 
normal strains to shear stresses. Thus, when Ai6 and A 26 are nonzero, and the laminate has a shear 
strain applied to it, normal stresses will result and vice-versa. These terms are analogous to the <2 1 6 
and 026 terms mentioned in the final part of section II. 

Equation (46) can be written as: 

b ,~\ £ [b#lm w )Mm)-£ ■ < 54 > 

L k=\ k= 1 z 

where t k is the thickness of the klh ply, and (h k +h k - 1)/2 is the distance from the geometric midplane 
to the center of the kth ply. This matrix is called the coupling stiffness matrix. From the constitutive 
equation (48), it can be seen that these terms relate bending strains (plate curvatures) with normal 
stresses and vice-versa. The B 16 and B 2 6 terms relate twisting strains to normal stresses and shear 

strains to bending stresses. 

If the laminate is symmetric, then the By terms will be the same for each mirrored ply above 
and below the midplane, with the exception of the sign of the (h k +h k - x )/2 term being negative if it is 
below the midplane (-z) and positive if it is above the midplane (+z). Thus, when summed, the result 

will be zero for all By. Now define: 

[ h k +h k-i] _ 


Part of equation (47) can be written as: 


[ h r h l-i) = [K + ^-i) ( h k- h k-l) + h\h k _ x -h k h\_ x 

= [K-^-i) 3 + 

= [(& k~h *_i) 3 + 3(* k -h k ^)(h k +h *-i) 2 -3(^-^L)] 



= t l +l2t h 


(55) 
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Therefore, equation (47) can be written as: 


n 

[5 

k= 1 








(56) 


It can be seen that the last term is the second moment of the kth ply with respect to the geometric 
midplane. D,y is called the bending stiffness matrix and relates the amount of plate curvatures with 
the bending moments. 


IV. NOMENCLATURE FOR DEFINING STACKING SEQUENCES 


A. Coordinate System 

The choice of coordinate system used for the laminate will determine its stacking sequence. 
For example, a unidirectional composite can be said to be made up of all 0° plies in the ^-direction, 
and all 90° plies in the y-direction. Alternatively, the same composite can be said to be made up of 
all 0° plies in the y-direction, and all 90° plies in the ^-direction. The composite can also be 
referenced by any other x-y coordinate system in the plane of the plate. The choice of coordinate 
system is totally arbitrary, but some general procedures are usually followed to make calculations or 
communication about the laminate easier for others to understand. 

A coordinate system is almost always chosen such that one of the axes runs in the direction 
of fibers of one of the plies of the laminate. This will make analysis much easier. The x-axis is 
usually chosen as the “longitudinal” axis, with the corresponding y-axis being the “transverse” 
direction. This is similar to what was defined earlier for laminae. The main load bearing fibers (if 
these are known) are usually called the 0° fibers, longitudinal fibers, or jc-direction fibers. The other 
ply orientations will then be defined with this coordinate system. An example of a typical coordinate 
system is given in figure 11. 


B. Nomenclature 

There is more than one way to denote the stacking sequence of laminates. However, once 
one method is learned, any other is easy to interpret, even though it may not be in the form that the 
user is accustomed to. 

Once the 0° fiber direction has been defined (and thus the x-axis), the plies that are not at 0° 
must be assigned an angle. To do this, start from the x-axis and rotate to the fiber direction of the 
ply being defined. Clockwise rotations are positive angles, and counterclockwise rotations are 
negative angles, although the reverse can also be used since only plane-stress is being examined for 
plates and the material is the same whether viewed from one surface or the other surface. Now that 
all plies have an angle associated with them, a method of presenting the stacking sequence follows. 

If the laminate is symmetric, then start with the angle of the outermost ply and write the ply 
angles, separated by a comma, until the midplane is reached. Enclose this string of angles in 
brackets or parentheses and subscript the brackets or parentheses with an “S” to denote 
“symmetric.” If the laminate is not symmetric, then proceed as above until the bottom ply is 
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Direction of primary loading 


1 


Laminate consisting of § plies 


+ 


4 


X 









Fiber orientation of the 8 plies 
as seen from the primary 
loading direction edge 



The stacking sequence from top 
to bottom is now 0°, +45°, -45°, 
90°, 90°, -45°, 445°, 0° 


y 


Figure 11. Coordinate system for a typical laminate. 
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reached. Subscript the brackets or parentheses with a “T” to denote “total” laminate. Referring to 
figure 10, this laminate is denoted as [0, +45, -45, 90]$ . This is more convenient than writing 
[0, +45, -45, 90, 90, -45, +45, 0] r . 

Further simplifications can be made when two or more plies of the same orientation are 
grouped together. The angle of these plies need only be written once with a subscripted number 
denoting the number of plies in the group. For example, [0, 90, 90, 90, 90, 0]r can be written as 
[0, 904 , 0]7’. Since this laminate is symmetric, further simplifications can result, and this laminate 
could be described by [0, 902]s- If a symmetric laminate consists of an odd number of plies, then the 
geometrical midplane of the laminate will lie at the midplane of the center ply. In this case, a bar is 
placed over the angle of this ply to denote that half of it resides in the top half of the laminate and the 
other half resides in the bottom half of the laminate. For example, a laminate with stacking sequence 
[0, 90, 90, 90, 0 ]t can be written as, [0, 90, 90] s . Any repeating units within the laminate can be 
placed in parentheses with a subscripted number representing the number of repeats. For example, a 
[0, 90, 0, 90, 0, 90, 0, 90]s laminate can be written as [(0, 90)4]$. If adjacent plies are of the same 
angle, but with different signs, then a plus-minus sign is usually placed in front of the angle of the 
plies. For example, a [0, +45, -45, 90, +30, -30]r laminate can be written as [0, ±45, 90 , ± 30] r- 
Some examples of stacking sequences and how they can be denoted are given in figure 12. 


445 


0 


0 

90 


0 


445 

-45 

[+45, 90, -45] s 

0 

[0 6 ] T 

-45 

-45 

0 

-45 

90 


0 


445 

445 


0 


0 


[0, ±45] s 


445 


-45 


90 

0 


90 

0 


90 

-45 

445 


[ ±45, 90, 0, 90] s 


0 


90 


0 


90 


0 


90 


0 


90 

0 


[(0,90) 2 , 0] s 


445 

-45 


0 

445 


-45 


0 

445 


-45 

0 


[(±45, 0)3 ) T 


Figure 12. Some laminate stacking sequences and their notation. 
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V. IN-PLANE ENGINEERING CONSTANTS FOR THE LAMINATE 


A. Orientation of the Laminate 

For a given stacking sequence of laminae whose engineering properties are known, it is 
possible to determine the in-plane engineering constants of the laminate from the Ay matrix for 
symmetric laminates, and the Ay, By, and Dy matrices for unsymmetric laminates. The choice of 
coordinates will determine the directions of the laminate engineering constants being evaluated. 
These directions are arbitrary, but are usually chosen as described in section IV. A. 


B. Symmetric Laminates 


Recall that for symmetric laminates, the By matrix consists of all elements being zero. This 
greatly simplifies finding the in-plane engineering constants of the laminate. To find the x-direction 
modulus, the value of the x- direction stress to the x-direction strain must be calculated. In equation 

form: 


E x = 



(57) 


where h is the thickness of the laminate. Since the Bys are zero, the constitutive equations are: 


> 


An Ai 2 Ai6 

V v 

= 

^12^22^26 



^16^26^66 


(58) 


Since a relationship between N x and £ x are being sought, when a load is applied in the x-direction, 

from equation (58): 


From equations (60) and (61): 


.. . o . o . o 

N x = A \\£ x +A n£ y +A \6r xy - 

()=■ A n £ x + A 22£y +A 2(,Yxy ’ 
^-A-i&x+AisEy+AteY^ ■ 



^26^16-A 12 A 66 
2 

^22^66“^ 26 J 


(59) 

(60) 
(61) 

(62) 
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and 


y = e 

>xy x 


_i4 16 ^ 26-^ 12^66“^ 26 ^ 16^ 

•* + 2 2 

^22-^66 “^26 ^66 J 


66 


( 63 ) 


Equations (62) and (63) can be substituted into equation (59) with the result: 


N x 

— =A n +A 


• 12 ' 


f A 26 A l6~ A l2 A 66^ 

2 

A 22 A 66~ A 26 


+ A 


16 


-A 16 A 26 A 12 A 66 -A 26 2 A 16 ^ 


^ 66 A 22 A 66 -A 26 A 


( 64 ) 


66 


Thus, E x can be calculated by dividing equation (64) by the thickness of the laminate, h which will 
give equation (57). 

The same procedure is followed to obtain E y . The constitutive equations are: 


Q~ A n£ x + A n £ +A l6 y , 


( 65 ) 


N y- A 12£ x +A 22 £ v +A 26Y xy ’ 


From equations (65) and (67): 


n < o . o . o 

0 ~ A \b e x +A 26 £ y +A 66^ ■ 


( 67 ) 
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x y 


A 16 ^ 26 ~^ 12^66 
AhA 6 6-Ai 6 2 ) 
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‘xy y 
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_ 26 | ^16^12^66"^ 16 ^26 
L 66 ^11^66 16 ^66 ) 


Equations (68) and (69) can be substituted into equation (66) with the result: 




0 - A 12l 


^ 16 ^ 26 ~^ 12^66 
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V ^11^66“^ 16 J 


+a 22 +a 26 
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~ A 26 ^16^12^66"^ 16^ 26 ^ 
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AllA 6 6 -A 16 A 66 j 
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( 70 ) 


Dividing equation (70) by the laminate thickness will give E y . 


29 



Gxy is found in the same manner. The constitutive equations are: 


From equations (71) and (72): 


« . 0 , 0 „ 0 
0 - A n e x +A 12 £ y +A , 


0 - A n£ x +A 22 £ y +A - 26Y xy > 


N xy- A \6 £ x +A 26 £ y +A 66Yxy • 


0 0 


A 12 ^ 26 -^ 16^2 2 
^ 11 ^ 22 ~^ 12 2 
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(71) 

(72) 

(73) 


(74) 


and 


r 


0 0 
£ y = ^ 


\ 


~ A 26 A 16 A 12 A 22~ A 12 A 26 

+ 5 2 

A 22 -A 1 1^22 “^12 ^22 


V 

Equations (74) and (75) can be used with equation(73) to obtain: 


) 


^xy _ A 26 ^ A \2 A \6 A 26 A 22~ A \2 A 26 ~ A 16 A 22 

0 ~ A<i6 ~ a + A A 2 . 2. 

7 xy Al2 A U A 22~ A 12 A 22 

Dividing equation (76) by the laminate thickness will give G^. 

To find Poisson’s ratio of the laminate, use equations (60) and (61) to obtain: 


0 - A \2 £ r +A 22 £ v + A 26 


A 16 0 A 26 0 


66 
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Rearranging to get: 


A 16 A 26 
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1 66 
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Using equations (65) and (67): 


Rearranging to get: 
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0 = A n £ x + A i 2 £ y +A u\ 
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A 66 x A 
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A 16^26 
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( 79 ) 


( 80 ) 


Example 1: 

Suppose a 4-ply laminate of AS4/3501-6 is laid up in a [0,+45]s stacking sequence. The ply 
properties are given below for this material: 

E x = 20,010,000 lb/in 2 , 

E 1 = 1,301,000 lb/in 2 , 

G 12 = 1,001,000 lb/in 2 , 

V 12 = 0.3 , 

v 2 i = 0.02 , 


Ply thickness = 0.005 in . 


This is all of the information needed to calculate the in-plane engineering constants. Calculating E x 
will demonstrate the general procedure applied to any of the engineering constants. From equation 
( 10 ): 

Q n = — = 2Q ; Q1Q { )QQ!b ^ — = 20,130,785 lb/in 2 


1 V 12 V 21 


1 -(0.3X0.02) 


222 “ 


Q 12 “ 


E, 


l ~ V 12 V 2l 


1,301,000 lb/in = 1>308;853 1Mn 2 ; 


1 -(0.3X0.02) 
v i 2 £ 2 (0.3)(1,301, 000 lb/in 2 ) 


l ~ v n v 2 i 


l-(0.3)(0.02) 


392,656 lb/in 


(El. 1) 


,. 2 


Q 66 = 1,001,000 lb/in . 


From equation (22), for the 0° plies: 
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Q° n = 20,130,785 lb/in 2 , 
Q° n = 392,656 lb/in 2 , 


0° 2 = 1,308,853 lb/in 2 , 
-0 

016-0 , 


-0 


2 26 “ ^ > 


0° 6 = 1,001,000 lb/in 2 . 


(El. 2) 


For the +45° plies: 


-445 


Q'n= (20,130,785 lb/in 2 )(0.707)%2( (392,656 lb/in z +2(l,001,000 lb/in 2 ) |(0.707) 2 (0.707) 


. 2, 


-+45 
2 12 


+ (1,308,853 lb/in 2 )(0.707) 4 = 6,557,237 lb/in 2 , 


(20, 1 30,785 lb/in 2 + 1 ,308,853 lb/in 2 -4( 1 ,00 1 ,000 lb/in 2 ) 
+392,656 lb/in 2 ((0.707) 4 +(0.707) 4 ) = 4,555,238 lb/in 2 , 


(0.707) 2 (0.707) 2 


-+45 


Q'22 = (20,130,785 lb/in")(0.707)‘+2^(392,656 lb/in 2 +2(l,001,000 lb/in 2 )j(0.707) 2 (0.707) 2 
+( 1 ,308,853 lb/in 2 )(0.707) 4 = 6,557,237 lb/in 2 , 


-+45 


Q 16 ~ = 20,130,785 lb/in 2 -392,656 lb/in 2 -2( 1,00 1,000 lb/in 2 ) 


. 2 


« 2 X 


(0.707) (0.707) 


392,656 lb/in"- 1,308,853 lb/in 2 +2(l, 00 1,000 lb/in 2 ) |(0.707)(0.707) J = 4,705,483 lb/in 2 , 


. 2 


-+45 


Q-26 = 20,130,785 lb/in 2 -392,656 lb/in 2 -2( 1,00 1,000 lb/in 2 ) |(0.707) 2i (0.707) 


■ 2 . 


2 


+ 392,656 lb/in"- 1,308,853 lb/in 2 +2(1, 00 1,000 lb/in 2 ) 


(0.707X0.707) 3 = 4,705,483 lb/in 2 , 


-+45 


20,130,785 lb/in 2 + 1,308,853 lb/in 2 -2(392,6561b/in 2 +l,001,000 lb/in 2 ) 
•1,001,000 lb/in 2 {(0.707) 4 +(0.707) 4 ] = 5,163,582 lb/in 2 . 


(0.707) 2 (0.707) 2 


(El. 3) 
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Now from equation (53), the Ay are obtained: 


. 2 


A n = [20,130,785+6,557 ,237+20,130, 785 +6,557,237] Ib/in (0.005 in) =266,880 T 


2 


A 12 = [392,656+4, 555, 238+392, 656+4, 555,238]lb/in (0.005 in) = 49,479 


m 


m 


A 2 2 = [1,308,853+6, 557,237+1,308, 853 +6,557,237] lb/in 2 (0.005 in) = 78,661 

A 16 = [0+4,705,483+0+4, 705, 483)lb/in 2 (0.005 in) = 47,055 — , 

in 

A 26 = [0+4,705,483+0+4, 705, 483]lb/in 2 (0.005 in) =47,055— , 

in 

A 66 = [1,001,000+5, 163,582+1, 001,000+5, 163 ,5 82] lb/in 2 (0.005 in) = 61,646 — 

in 


(El. 4) 


From equations (64) and (57): 
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E = 


lb 


in 
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(49479) 


(47,055)(47,055)-(49,479)(6 1,646) 
(78,66 1)(6 1 ,646)-(47,055)" 
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in 
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(47,055) 


- (47,055) 
(61,646) 
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(47,055)(49,479)(61,646)-(47,055) 2 (47,055) 

(78,66l)(61,646) 2 -(47,055) 2 (61,646) 


lb 

in 


0.02 in 


E = 1 1,333,000 lb/in 2 


(El. 5) 


C. Nonsymmetric Laminates 

Since for a nonsymmetric laminate there are nonzero By terms, the calculation of in-plane 
engineering constants become more involved; however, the same basic procedure is followed as for 
symmetric laminates. Since there are now six equations instead of three, it will be much easier to 
use matrix notation. The constitutive equations are given by equation (48): 



N x 


An A 12 A 16 

B it B 12 B i6 


£° 

n 

N y 


A 12 A 22 A 26 

B 12 B 22 B 26 
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A 16 A 26 A 66 

B 16 B 26 B 66 


yO 
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M x 
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B n B 12 B i6 
B 12 B 22 B 26 
B 16 B 26 B 66 

£*11 *>12 £>16 
£>12 £*22 £>26 
£>16 £>26 £>66 
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To find E x , only the x-direction in-plane load is applied and a relationship between N x and e® is 
sought. The constitutive equation now becomes: 
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Using Cramer’s rale to solve for e\ : 
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B 16 B 26 B ^6 
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£>26 
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£>26 
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( 81 ) 


( 82 ) 
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Determinants for two 6x6 matrices must be found. 


At this point, it should be clear why a computer (or calculator) program makes calculations 
much easier. To write out the solution for equation (82) would require 6! = 720 terms for the 
numerator and the denominator. Cofactor expansion can be used in the numerator for some 
simplification: 
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From equation (57), E x can be found by evaluating: 


N x /h 
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Appendix A shows a Fortran program to calculate the determinant of a 6x6 and 5x5 matrix. 
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Ey can be found in a similar manner. The denominator will be different since equation (81) is 

being solved for e ® . 


Gxy will be given by: 
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Poisson’s ratio will be determined as it was for symmetric laminates. For y*y, the contraction in the 
y-direction upon an applied stress in the jc-direction must be obtained. This is given by: 
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Thus, Vxy will be given by: 
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Vy* is given by: 
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Example 2: 

Given a 2-ply laminate of stacking sequence [0,+45]y, E x will be calculated. From example 1, 
the Ay’s are known. The Bi/s and Dy’s must be calculated. 


To obtain the By, use equation (54) with the already determined values in equations (El. 2) 

and (ELS). 


B n = 20,130,755 lb/in 2 (0.005 in)(0.0025 in) 
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B 16 = 0+4,705,483 lb/in 2 (0.005 in)(-0.0025 in) = -59 lb , 

B 26 = 0+4,705,483 lb/in 2 (0.005 in)(-0.0025 in) = -59 lb , 

B 66 = 1,001,000 lb/in 2 (0.005 in)(0.0025 in) 

+ 5,163,582 lb/in 2 (0.005 in)(-0.0025 in) = -52 lb . 


Use equation (56) to find the D/s: 
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Now use equation (84) and the Fortran programs in appendix A to obtain: 
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D. Summary 


The equations for the in-plane engineering constants of a symmetric laminate are 
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For unsymmetric laminates, the engineering constants are given by: 
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VI. ENVIRONMENTAL EFFECTS 


A. Importance 

Like all engineering materials, laminates will contain residual stresses from processing, or 
will experience strains due to thermal effects or moisture absorption. However, because of the 
anisotropic nature of composite laminates, these effects become much more important. It is well 
known that most laminates are processed or “cured” at elevated temperatures. It is at this elevated 
temperature that the molecular structure of the material is set. Upon cooling, the laminate will 
experience many internal stresses as each ply contracts a different amount in different directions. 

For unsymmetric laminates, this can be seen by warpage of a flat plate as it is removed from its 
molding platens. In some cases, this warpage is desirable to achieve a natural “twist” in the 
material (such as for helicopter rotors). In general, the more anisotropic the laminate, the more 
important the residual thermal and moisture absorption stresses become. The thermal effects will be 
quite important in the next section on ply stresses. 


B. Coefficients of Thermal Expansion 

Just like any other material, fiber/resin systems will experience a change in strain with a 
change in temperature. These strains are defined by the coefficients of thermal expansion. These 
values are material constants in each principal material direction. Thus, two constants will describe 
the thermal expansion coefficients for any lamina. These are defined as: 

cq - coefficient of thermal expansion in the fiber direction 


and 


a 2 - coefficient of thermal expansion in the direction perpendicular to the fibers. 


43 



These values have dimensions of inch/inch/°F or inch/inch/°C, depending on the temperature scale 

that is being used. 

Just as with the mechanical strains, the thermal strains must be transformed into the 

laminate coordinate system: 
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Note that a coefficient of thermal expansion for shear is formulated if the lamina being examined does 
not have its material axes as principal material axes (i.e., 0° or 90°). The amount of thermal 

sixain induced in each lamina is given by: 

e x = a x AT , 


a y AT , 


T = lk 

£ xy 9 


(XxyAT , 


(91) 


where the superscript T denotes “thermal,” and AT denotes the change in temperature from cure to 

operating temperature. 

These thermal strains are now treated just like the mechanical strains considered earlier. 
Therefore, from equation (43), it can be seen that when these thermal strains are combined for each 
layer of the laminate, thermal stress resultants are present: 
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From equation (44), it can be seen that thermal moment resultants also develop: 
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These thermal stress and moment resultants can be added to the mechanical stress and moment 
resultants to arrive at the total stress and moment resultants: 
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where the superscript Tot denotes “total”. 

Typical values of the thermal expansion coefficients for carbon/epoxy is a 1 = -0.07 2x1 0- 6 
in/in/°F and a2 = 32.4x1 0~ 6 in/in/°F. Note that in the fiber direction, the material actually contracts 
upon heating (or expands upon cooling). Also, note the large differences in the two coefficients. For 
this reason, composite materials have been used in many applications where a certain thermal 
expansion is desired, since plies can be combined to give a wide range of values (including near zero 
which is very convenient for optical benches). 


C. Moisture Effects 

Swelling of a composite material due to moisture absorption is handled in the exact same 
manner as expansion due to temperature differences. The moisture swelling coefficients in each 
principal material direction must be known. They are designated as pi, along the fiber direction, and 
p2, perpendicular to the fiber direction. The strain due to moisture absorption is given by: 

t 

M tl A 

=P X A m , 

Sy=P y Am , (95) 

M a a 

Y - p Am , 

'xy i-’xy ’ 

where the /Vs are the transformed moisture expansion coefficients, Am is the moisture concentration 
in weight moisture/weight material, and the superscript M denotes “moisture.” 

Stress and moment resultants can be determined due to the effects of moisture absorption. 

Typical values of moisture expansion coefficients for carbon/epoxy are, pi = 0.01 in/in/g/g and 
p2 = 0.35 in/in/g/g. The moisture concentration is usually a very low number under normal operating 
circumstances -0.0005 g/g. However, in humid environments this number may be much higher. 
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VII. STRESSES AND STRAINS WITHIN LAMINAE OF A SYMMETRIC LAMINATE 


A. Strains Within the Laminae 

So far, in this paper, methods have been shown to describe the stress-strain behavior of a 
laminate as a whole. In many instances, however, the stresses and/or strains in each ply may need 
to be known. To do this, the constitutive equation (52) is used to find the total midplane strains and 
curvatures. From these values, the total strain in each ply can be calculated from equation (30). 
These strains can then be transformed into strains in the principal material directions by the first 
equation in (18). As in finding the engineering constants of a laminate, the most involved step in 
performing this analysis is the manipulation of matrices. The A', B', and D' matrices must be found 
from the A, B, and D matrices. 


Example 3: 

Using the symmetric laminate in example 1, [0,4-45]$, find the strains in each ply in each ply’s 
principal material direction given a tensile load of 1,000 lb on a 2-in wide specimen in the ^-direction. 
The x-direction is parallel to the 0° fibers. Assume no environmental effects. 

The specimen is 2-in wide, therefore, the stress resultant is 
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From example 1, the Ay are known. Since the laminate is symmetric, the By are zero. The £>y are 
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Using equations (51) and (53), the A' and D' matrices can be found. The B' and C matrices will be 

zero since the B matrix is zero. 
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Therefore, the constitutive equation is: 
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These are the midplane strains and curvatures. Using equation (30), it can be seen that since the 
midplane curvatures are zero, then the strains in each ply will be equal to the midplane sixains. In 
other words, the strains are constant through the thickness in any given direction since there is no 
bending and therefore z does not enter into the calculation. 

To transform these strains into the principal material directions for the +45° ply: 
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Note: Recall that the engineering shear strain must be put into tensorial shear strain, % = 1/2 Yxy, 

before transformation. 


B. Stresses Within the Laminae 


The stresses in a lamina can be determined once the strains are known by applying equation 
(23) to each ply. If thermal or moisture strains are present, these strains must be subtracted from 
the total strain since they are not caused by an external force. 


Example 3a: 

Find the stresses in each ply of the laminate in example 3. 


<7 

X 


20,130,785 

392,656 

0 


0.00221 


<*x 


44,214 



= 

392,656 

1,308,853 

0 

lb/in 2 

-0.0007 


°y 

= 

48 

lb/in 2 

y 

X 

xy 

0 

0 

0 

1,001,000 


-0.00115 

0 


0 

-1,151 



(E3a.l) 




6.557.237 4,555,238 4,705,483 

4.555.238 6,557,237 4,705,483 

lb/in 2 

0.00221 

-0.0007 


°x 

°y 


5,892 

66 

T xy 

+45 

4,705,483 4,705,483 5,163,582 


-0.00115 

+45 


+45 

1,167 


lb/in 2 


(E3a.2) 


Transforming the +45° ply into its principal material directions: 


<*\ 


(j 2 

— 

T 12 

-f45 

[ 


cos 45° 
sin 2 45° 


sin 2 45° 


cos 45° 
sin 45° cos 45° sin 45° cos 45° 


2sin 45° cos 45° 
-2sin 45° cos 45° 

cos 2 45° - sin 2 45° 


5892 

66 

1167 


lb/in 


o-f 5 = 4,146 lb/in 2 ; 


crj 45 = 1,812 lb/in 2 ; = -2,9 1 3 lb/in 2 


(E3a.3) 
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Example 4: 


The importance of thermal effects can be found by working problem 3, taking into account that 
the laminate was cured at 176.7 °C and tested at 21.1 °C. This implies a AT of -155.6 °C. The 
coefficients of thermal expansion for AS4/3501-6 are given below: 

a, = -0.04x1 0~ 6 ; cc, = 18xl0 -6 . (E4.1) 

1 in °C 2 in °C 

Using equation (90), for the 0° plies the coefficients are: 

a = -0.04x1 0~ 6 ; a= 18xl0 -6 -^- , (E4.2) 

x in °C y in °C 


and for the +45° plies: 


+45 o no 1A~ 6 in 445 O AO 1A~ 6 m 445 AAA 1 A“ 6 111 

a =8.98x10 ; a =8.98x10 ; a =-9.02x10 

* in°C y in °C ^ in°C 

From equation (91), the thermal strains in the 0° and the 45° plies can be found: 


(E4.3; 


or m o t or 

£ =0.0000062- ; e =-0.0028 ; y =0 

jc in y ,x y 


(E4.4) 


e +45T = -0.0014— ; s 45T = -0.0014 ; y +45T = 0.0028 (E4.5) 

x in y ^ 

Using equations (92) and (93), these can be converted into thermal stress and moment resultants: 



T 

N x 


-33.57 


n t 

y 

= 

-60.42 


T 

N 


12.83 

The thermal moment resultants will be zero for the same 


in 

lb 


(E4.6) 


symmetric. Thus: 
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^ 11 ^ 12^16 
/ / r 

AAA 
/i 12 /1 22^26 


0 0 0 
0 0 0 


A 16 A 26 A <,<, 0 0 0 


o 0 0 D n D n D l6 


0 0 0 ^ 12 ^ 22^26 

0 0 ° D 'l6 D 26 D 66 


. 2 

.-6 in 


e = 4.41x10 — — 466.43- -1.4x10 — 
x lb in I lb 


. 2 
-6 in 


-60.42 -] -2.3X10 -6 — ( 12.83 — ] = 0.0021 1 
in | lb 1 in / 


g° = -1.4xl0~ 6 — (466.43-] +23.8xl0~ 6 —(-60.42 -] -17.1xl0 -6 — ( 12.83 - ] = -0.00231 
y lb in J lb ^ in J lb ^ in J 

y = -2.3xl0~ 6 — 466.43 -] -17. lxlO -6 — (-60.42 -] +31.05X10 -6 — ( 12.83 -] = 0.00036 
lb in J lb ^ in] lb { m J 

K x = 0 ; K y = 0 ; K^ = 0 


Transforming these into principal material directions for the 45° plies gives: 


£ l 


e 2 

= 

r , 2 


2 

+45 


cos 45° 
sin 2 45° 


sin 45° 
cos 2 45° 


2sin 45° cos 45 c 


0.00211 


sin 2 45° cos 45° -2sin 45° cos 45° -0.00231 

- sin 45° cos 45° sin 45° cos 45° (cos 2 45° - sin 2 45°] 0-00036 


0.0001 ; e 2 = -0.0003 


= - 0.00442 


The lamina stresses are: 
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°x 


20,130,785 

392,656 

0 


0.00211-0.000006 

G y 

= 

392,656 

1,308,853 

0 

lb/in 2 

-0.0023 1-(-0.0028) 

T 


0 

0 

1,001,000 


0.00036-0 

T *y 

0 

- 


■ 


■ 


°x 


42,543 

°y 

= 

1,467 

T 


360 

xy 

0 



For the 0° ply, the x-direction, stress is altered slightly, but the y-direction and shear stresses are 
quite different due to the thermal effects. 


<*x 


6,557,237 4,555,238 4,705,483 


0.0021 1 -(-0.0014) 

°y 

= 

4,555,238 6,557,237 4,705,483 

lb/in 2 

-0.0023 1-(-0.00 14) 

T 


4,705,483 4,705,483 5,163,582 


0.00036-0.0028 

Z xy 

+45 

- 


- 


J+45 


°x 


7,389 

<fy 

= 

-1,460 

T 


-366 

xy 

+45 



lb/in 


(E4.10) 


Transforming the +45° ply stresses into principal material direction stresses: 




*— H 

b 


^2 

= 

CN 

+45 


cos 45° 
sin 2 45° 


sin 2 45° 
cos 2 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 


sin 45° cos 45° sin 45° cos 45° cos 2 45° - sin"45° 


dl 45 = 2,600 lb/in 2 


c-! 45 = 3,330 lb/in 2 


7,389 

-1,460 

-365 


lb/in 2 


T ^ 5 = -4,425 lb/in 2 


(E4.ll) 


For the + 45° plies, the stresses vary substantially in all directions due to thermal effects. 


Example 5: 

Perform the same exercise as above taking into account moisture uptake. The amount of 
moisture absorbed by the composite is 0.007 g/g, and the moisture swelling coefficients are pi = 0.01 
in/in/g/g and = 0.35 in/in/g/g. 
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The moisture expansion coefficients need to be transformed for the +45° plies; 


P x 

Py 


A 




+45 


cos 45° 
sin 2 45° 


sin 2 45° 


-2 sin 45° cos 45° 


cos 45° 2 sin 45° cos 45 c 


sin 45° cos 45° - sin 45° cos 45° cos 45° - sin 45 c 


0.01 

0.35 

0 


<iT 45 = 0.18 in/in/g/g ; /f 45 = 0.18 in/in/g/g ; 0^ 5 = -0. 17 in/in/g/g 
The strains due to moisture uptake are given by equation (95): 


(E5.1) 


s° M = 0.00007— ; e° M = 0.00245— ; y M = 0 
* in y in ** 


(E5.2) 


e +45M = 0.00126 — . £ +45M = 0.00126— ; y U5M = -0.00238 — 
x in y in ^ in 


The stress and moment resultants due to moisture swelling can now be calculated using equations 

(92) and (93): 

VI 




51.7 

60.4 

-4.3 


in 

lb 


(E5.2) 


The moisture moment resultants will be zero since the laminate is symmetric. 

The new total stress resultant is the sum of the mechanical, thermal, and moisture stress 

resultants: 


N 

N 

N 


Tot 

x 

Tot 

y 

Tot 


500 

0 

0 



-33.6 

-60.4 

12.8 


IS. 

lb 


51.7 

60.4 

-4.3 


HI - 
lb “ 


518.1 

0 

8.5 


in 

lb • 


(E5.4) 


Using the constitutive equation, the midplane strains can be calculated as they were in equations 

(E3.4) and (E4.7). 
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y 

o 


7r 


xy 


K x 


K 


^ 11 ^ 12^16 0 0 0 
0 0 0 
0 0 0 


4 4 4 

12 22 26 


^ 16 ^ 26^66 


0 0 0 


D n D n D i6 


. 2 

o . . 1 . n — 6 in 
s =4.41x10 — 
* lb 

• 2 f 

0 -t a i r\~6 

e =-1.4x10 — 
y lb 




( ]h \ 

518.1 — 

V in 


0 0 0 ^ 12 ^ 22^26 
0 ° 0 


N. 


Jot 


N, 


X 

Tot 


N. 


y 

Tot 


xy 

0 

0 

0 


-1.4x10 


. 2 
-6 m 

lb 


518.! — 
in 


+23.8x10 


. 2 
-6 in 

lb 


f lb' 

. 2 
in. 

0 c lh \ 


-2.3x10 — 

8.5— = 

l in . 

lb 

k “i 

lb\ 

. 2 
— 6 in 

f lb^ 

0 — 

-17.1x10 — 

8.5 — 

in/ 

lb 

l in J 


= -0.00087 


0 O -2 1 A- 6 W 2 

y =-2.3x10 — 
Ix y lb 


518.1 —) -17.1x10 6 — 
in / lb 


^ Ibl „„ __ . _ -6 in 2 f _ lb 
0— +31.05x10 — 8.5 — 
in lb l in 


= -0.00093 


K x = 0 ; K v = 0 ; K xy = 0 


Transforming these strains into principal material direction strains for the +45° plies; 


£ i 


cos 2 45° 

sin 2 45° 

£ 2 

= 

sin 2 45° 

cos 2 45° 

Zh 


- sin 45° cos 45° 

sin 45° cos 45° 

2 

+45 

■ 



+45 


+45 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

cos 2 45° - sin 2 45° 


+45 


=0.00024 ; =0.00117 ; 7 ^ =-0.00314 

The stresses in the 0° plies are given by: 


0.00227 

-0.00087 

-0.00093 


. (E5.5) 


(E5.6) 
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f ] 

°x 


20,130,785 

392,656 

0 


0.00227-0.000006-0.00007 

CT V 

— 

392,656 

1,308,853 

0 

lb/in 2 

-0.00087-(-0.0028)-0.00245 

y 

X 

L X:y 

0 

0 

0 

1,001,000 


-0.00093-0-0 




43,963 

a 

= 

181 

y 

X 


-931 

xy 

0 



(E5.7) 


The stresses in the +45° plies are given by: 


K 


6,557,237 4,555,238 4,705,483 


0.00227-(-0.00 1 4)-0.00 1 26 



~ 

4,555,238 6,557,237 4,705,483 

lb/in 2 

-0.00087-(-0.0014)-0.00 1 26 


y 

X 

xy 

+45 

4,705,483 4,705,483 5,163,582 


-0.00093-0.0028-(-0.00238) 

+45 




6,125 

a 

= 

-161 

y 



T 


934 

T *y 

+45 

L 


(E5.8) 


Transforming these stresses into principal material directions for the +45° plies: 


u 2 

r !2 


J+45 


cos 2 45° 
sin 2 45° 


sin 2 45° 
cos 2 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

\ 


sin 45° cos 45° sin 45° cos 45° |cos 2 45° - sin 2 45° j 


6,125 

-161 

934 


lb/in 2 


+4S 9 

cr = 3,916 lb/in 


(E5.9) 


The swelling strain of moisture uptake tends to diminish the contraction due to processing thermal 

differences. 


Example 6: 

Find the stresses and strains in each lamina of the laminate used in the previous examples if 
the laminate is subjected to a bending moment in the 0° fiber direction of 5 in-lb/in. 
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The constitutive equation is 


0 


/ / / 
^11^12^16 

0 0 0 

0 


t / / 




A 12 A 22 A 26 

0 0 0 

0 


/ / / 

0 0 0 



A 16 A 26 A 66 

K x 


0 0 0 

^11 ^12^16 

Ky 


0 0 0 

D n D 22 D 26 

K 

/v xy 


0 0 0 

D 16 D 26 D 66 


0 

0 

0 

M x 

0 

0 


K =D. .M = 0.0836 7—77 
11 m-lb 


in-lb 


in 


= 0.418 


m 


K y = D 12 M x = - 0.033 7 


in-lb 


K xy = D.M x = ~0M95—r 
^ 16 x in-lb 


f 5 in ~ lb ^ 
in , 

in-lb 


= -0.165 r 


in 


) 

£ ° = 0 ; e° = 0 ; y =0 ' 

x y 'xy 


in 


= -0.0975 - 
in 


(E6. 1 ) 


Using equation (30), it can be seen that the strains in each ply will vary across the thickness 
of the ply. Thus, the distance the ply is away from the geometric midplane must be taken into 
account, as well as the direction of the fibers as in the previous examples. For clarity, assign the four 
plies of the laminate numbers as shown in the following: 


Z 



Since this laminate is symmetric, the geometric midplane is also the neutral plane of the 
plate. Thus, the strains will be zero at the bottom of ply No. 2 and at the top of ply No. 3. The 
following are calculations for other planes within the plate. 
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The top of ply No. 2: 


- 

Top 

1 



e , 

£ 

= 0.005 in 

0.418 

-0.165 


0.00209 

-0.00083 

y 

v 


-0.0975 

in 

-0.00049 

fxy 

+45 



L 


(E6.2) 


The bottom of ply No. 1 is at the same geometric location as the top of ply No. 2 (i.e., at z- 0.005 
in). Thus, the strains will be the same. 


The top of ply No. 1: 


£ 

Top 

= 0.01 in 

0.418 

-0.165 


0.00418 

-0.00165 

y 

0 

-0.0975 

in 

-0.000975 


The bottom of ply No. 3: 


£ 

Bot 

-0.005 in 

0.418 

-0.165 


-0.00209 

0.00083 

y 

Yxy 

+45 


-0.0975 

in 

0.00049 


(E6.3) 


(E6.4) 


The top of ply No. 4 is at the same location as the bottom of ply No. 3 and, thus, will have the 

above strains. 


The bottom of ply No 4: 


£ 

Bot 

-0.01 in 

0.418 

-0.165 


-0.00418 

0.00165 


0 


-0.0975 

in 

0.000975 


(E6.5) 



Transforming these strains into strains in the principal material directions for the +45 ° plies: 
Top of ply No. 2: 






£ 2 

= 

Tl2 


2 

+45 


cos 2 45° 
sin 2 45° 

- sin 45° cos 45° 


sin 2 45° 
cos 2 45° 

sin 45° cos 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

f ' 

cos 2 45° - sin 2 45° 


0.00209 

-0.00083 

-0.00049 

2 


4 - 4-5 4-45 

=0.00039 ; £2 = 0 - 00088 


7^2 = -0.00292 


Bottom of ply No. 3: 





e i 


cos 2 45° 

£ 2 

= 

sin 2 45° 

Yn 


- sin 45° cos 45° 

2 

+45 

• 


sin 2 45° 
cos 2 45° 

sin 45° cos 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

|cos 2 45° - sin 2 45° j 


e ^ 5 = -0.00039 ; e £ 5 = -0.00088 ; 7^ = 0.00292 


-0.00209 

0.00083 

0.00049 

2 


The stresses within the laminate is calculated from the strain using equation (23): 
Top of ply No. 1: 


1 

Top 

r 


1 


0.00418 

°X 


20,130,785 

392,656 

0 


°y 

= 

392,656 

1,308,853 

0 

lb/in 2 

-0.00165 

T 


0 

0 

1,001,000 



xy . 

0 



• 


-0.000975 


• 

Top 




83,500 

°y 

= 

-518 

T 


-976 

xy 

0 



(E6.6) 


(E6.7) 


(E6.8) 
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Bottom of ply No. 1: 


■ 

Bot. 

<T 


X 


<7 

= 

y 


X 



0 


20,130,785 

392,656 

0 


0.00209 

392,656 

1,308,853 

0 

lb/in 2 

-0.00083 

0 

0 

1,001,000 


-0.00049 


• 

Bot 


°x 


41,747 

a 

= 

-266 

y 



T 


-490 

T xy 

0 

L 


lb/in 2 


Top of ply No. 2: 


<*x 

<y v 

Top 

6.557.237 4,555,238 4,705,483 

4.555.238 6,557,237 4,705,483 

lb/in 2 

0.00209 

-0.00083 

y 

> 

+45 

4,705,483 4,705,483 5,163,582 


-0.00049 

. 


+45 


- 

Top 

1 

°x 


7,618 

°y 

= 

1,772 

X 


3,399 

xy 

+45 



(E6.9) 


(E6.10) 


Since the laminate is symmetric, the stresses on the bottom half of the laminate will mirror 
those on the top, only with a change of tensile stresses to compressive and compressive to tensile 

(i.e., a change in sign). 

Transforming the stresses into the principal material directions for the +45° ply: 


Top of ply No. 2: 


(J i 

Top 

^2 

= 

T 12 

+45 


cos 2 45° 
sin 2 45° 

sin 45° cos 45° 


sin 2 45° 
cos 2 45° 

sin 45° cos 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

|cos 2 45° - sin 2 45°j 


7,618 

1,772 

3,399 


lb/in 2 



Top 

8,094 

*2 

= 

1,296 

T 12 

+45 

-2,923 


lb/in 2 


(E6.ll) 
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The stresses in the bottom of ply No. 3 will be the mirror of the stresses in the top of ply No. 2, thus, 
for the bottom of ply No. 3: 


*1 

Bot 

-8,094 

*2 

T 12 

+45 

-1,296 

2,923 


lb/in 2 


(E6. 12) 


All of the information in example No. 6 can be tabulated: 


Strains 


Plane 

z (in) 

£ x 

£y 

Yxy 

£ i 

e 2 

ft 2 

Ply 1 (Top) 
Ply 1 (Bot.) 

+0.01 

+0.005 

0.00418 

0.00209 

-0.00165 

-0.00083 

-0.00098 

-0.00049 

0.00418 

0.00209 

-0.00165 

-0.00083 


Ply 2 (Top) 
Ply 2 (Bot.) 

+0.005 

0 

0.00209 

0 

-0.00083 

0 

-0.00049 

0 

0.00039 

0 

0.Q0088 

0 

H9 

Ply 3 (Top) 
Ply 3 (Bot.) 

0 

-0.005 

0 

-0.00209 

0 

0.00083 

0 

0.00049 

0 

-0.00039 

0 

-0.00088 

0 

0.00292 

Ply 4 (Top) 
Ply 4 (Bot.) 

-0.005 

-0.01 

-0.00209 

-0.00418 

0.00083 

0.00165 

0.00049 

0.00098 

-0.00209 

-0.00418 


0.00049 

0.00098 


Stresses 


Plane 

z (in) 

o x 

°y 

^xy 

<y\ 

<r 2 

?12 

Ply 1 (Top) 

+0.01 

83,500 

-518 

-976 


-518 

-976 

Ply 1 (Bot.) 

+0.005 

41,747 

-266 

-490 

41,747 

-266 


Ply 2 (Top) 

+0.005 

7,618 

1,772 

3,399 


1,296 

-2 923 

Ply 2 (Bot.) 

0 

0 

0 

0 


0 

WM.'? s 

Ply 3 (Top) 

0 

0 

0 

0 


0 

0 

Ply 3 (Bot.) 

-0.005 

-7,618 

-1,772 

-3,399 


-1,296 

2,923 

Ply 4 (Top) 

-0.005 

-41,747 

266 

490 

-41,747 

266 

490 

Ply 4 (Bot.) 

-0.01 

-83,500 

518 

976 


518 

976 


Since the stresses and strains vary linearly through the thickness of each ply, a diagram of the 
stresses and strains through the laminate can be made. A few examples follow. 
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VIII. STRESSES AND STRAINS WITHIN LAMINAE OF AN 
UNSYMMETRIC LAMINATE 


A. Difference From Symmetric Laminates 

If the laminate being examined is not symmetric, then additional complexities arise in the 
material’s behavior from that examined in the previous section. The neutral plane of the plate will not 
coincide with the geometrical midplane of the plate as it did for symmetric plates. Thus, there will be 
less simplifications, such as knowing the bottom half of a symmetric laminate will consist of a 
negative mirror image of stresses and strains from the top half due to bending moments. Also, 
environmental effects may become critical since the plate can warp due to these effects. Coupling 
now becomes extremely important and can cause very unique mechanical behavior characteristic of 
anisotropic plates 


B. Example (0/+45)j Laminate 

Example 7: 

Consider a 2-ply laminate of AS4/3501-6 with a (0/+45)7- stacking sequence. The lamina 
properties were given in example 1 as: 


Ei = 20,010,000 lb/in 2 , 
E 2 = 1,301,000 lb/in 2 , 
G 12 = 1,001,000 lb/in 2 , 
V 12 = 0.3 , 
v 2 i =0.02 , 

Ply thickness = 0.005 in , 


and the resulting values for Q- for the 0° and the +45° plies were given as: 


Q° n = 20,130,785 lb/in 2 , 

— h45 2 

Q n = 6,557,237 lb/in , 

Q° n = 392,656 lb/in 2 , 

— h45 2 

<2 12 = 4,555,238 lb/in , 

Q° 12 = 1, 308,853 lb/in 2 , 

Q 22 = 6,551‘. 237 lb/in 2 , 


Q 16 = 4,705,483 lb/in 2 , 

226 = ° • 

— +45 2 

Q 26 =4,705,483 lb/in , 

Q° 66 = 1,001,000 lb/in 2 , 

<2 66 =5,163,582 lb/in 2 . 
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From equation (53), the Ay are obtained: 


A n = [20,130,785+6,557,237] lb/in 2 (0.005 in) = 133,440— , 

in 

A n = [392,656+4,555,238] lb/in 2 (0.005 in) = 24,739 — , 

in 

A 22 = [1,308,853+6,557,237] lb/in 2 (0.005 in) =39,330 — , 

in 

A 16 = [0+4,705,483] lb/in 2 (0.005 in) = 23,527 — , 

in 

A 26 = [0+4,705,483] lb/in 2 (0.005 in) =23,527 — , 

in 

A 66 = [1,001,000+5,163,582] lb/in 2 (0.005 in) = 30,823 — . 

in 

From equation (54), the By can be determined: 

B n = [20,130,785-6,557,237] lb/in 2 (0.005 in) (0.0025 in) = 170 lb , 
5 12 = [392,656-4,555,238] lb/in 2 (0.005 in) (0.0025 in) = -52 lb , 
B 22 = [1,308,853-6,557,237] lb/in 2 (0.005 in) (0.0025 in) = -66 Jb , 
B l6 = [0-4,705,483] lb/in 2 (0.005 in) (0.0025 in) = -59 lb , 
B 26 = [0-4,705,483] lb/in 2 (0.005 in) (0.0025 in) = -59 lb , 
B 66 = [1,001,000-5,163,582] lb/in 2 (0.005 in)(0.0025 in) = -52 lb . 


From equation (56), the Dy can be determined: 


D„= 20,130,785 + 6,557,237 


lb/in 


(o.005 in.)' 
12 


+ { 0.005 in.)(0.0025 in.)' 


= 1.11 in-lb , 


J 


D 


12 ’ 


392,656 + 4,555,238 


lb/in 2 


(o.005 in.) 3 
12 


+ 


(o.005 in.)(o.0025 in.)' 


(E7.1) 


(E7.2) 


= 0.206 in-lb , 



= 0.328 in-lb , 


D 22 = 


1,308,853 + 6,557,237 


lb/in 2 


(o.005 in.) ( y y 

± ■^- L - + (0.005 in.)(0.0025 in.) 


£>16 = 


0 + 4,705,483 


lb/in 2 


(o.005in.) ( y y 

± + (0.005 in.)(0.0025 in.) 


= 0.196 in-lb 


J 


!( V 

r -1 2 ( 0.005 in.) / y a 

D 26 = [o + 4,705,483J lb/in 4 ___2_ + ( 0 .005 in.)(0.0025 in.) 


£>66 ~ 


1,001,000 + 5,163,582 


lb/in 


( 0.005 in.) t y 

1 ’ + (0.005 in.)(0.0025 in.) 


0.196 in-lb 
A 


V 


(E7.3) 


= 0.257 in-lb . 


J 


Find the stresses and strains in each of the two plies given a stress resultant in the x-direction of 
250 lb/in. The primed matrices must be obtained from equations (51) and (53): 



8.80 -2.80 -4.60 
-2.80 47.7 -34.3 
-4.60 -34.3 62.1 


in— 6 m 

x 10 — 
lb 



-1.91 0.0014 0.115 
0.933 0.979 0.866 
02.66 1.16 0.934 


x 10 3 in , 



1.91 -0.933 -2.66 

-0.0014 -0.979 -1.16 
-0.115 -0.866 -0.934 


x 10 3 in , 



0.580 0.087 0.115 
0.087 0.195 0.078 
0.115 0.078 0.151 


in-lb . 


Therefore, 




17.0 

-3.07 

-11.6 

-3.07 

54.6 

-25.9 

-11.6 

-25.9 

77.7 


x 10 


-6 in 
lb 



-4.07 

0.353 

3.68 

0.368 

3.38 

3.72 

3.75 

3.72 

1.43 




-4.07 0.368 
0.353 3.38 
3.68 3.72 


3.75 

3.72 

1.43 




2.05 

-0.366 

-1.37 

-0.366 

6.53 

-3.09 

-1.37 

-3.09 

9.27 


Using the constitutive equation (52), the midplane strains and plate curvatures can be found. 

( 


e=a'n x = 17.0x10 6 — 
* 11 * lb 


250 — | = 0.00425 
in 


e° = a'n x = -3.07x10 6 — 

y 12 * lb 


f 


250 — I = -0.00077 
in 


r° =a’n x = -11.6x10 6 — 

'xy 16 x p, 


f 


250 — | =-0.0029 
in 


Kx=C„N x = - 4.07x10 3 - 
x 11 lb 


V 


1 1U ) 

250 — 
in 


= - 1.02 


/ 


m 


(E7.4) 


K=C. 0 N= 0.353x10 - 
y u lb 


250-1 = 0.088- 


V 


in , 


in 


K xv =C u .N= 3.68x10 - 
*y lo •* ||-) 


/ 

250 — 


in 


= 0.92 


in 


66 



The strains throughout the laminate can now be calculated from equation (30). Use the following 
diagram for clarity; 



At the top of ply No. 1: 


£ 

Top 

0.00425 

-0.00077 

+ 0.005 in 

-1.02 

0.088 


-0.00085 

-0.00033 

y 

0 

-0.0029 


0.92 

in 

0.0017 


At the bottom of ply No. 1 : 


e 

Bot 

0.00425 

-0.00077 

+ 0 in 

-1.02 

0.088 


0.00425 

-0.00077 

^ & 

0 

-0.0029 


0.92 

in 

-0.0029 


At the top of ply No. 2: 


£ 

Top 

0.00425 

-0.00077 

+ 0in 

-1.02 

0.088 


0.00425 

-0.00077 


+45 

-0.0029 


0.92 

in 

-0.0029 


At the bottom of ply No. 2: 


£.. 

Bot 

0.00425 

-0.00077 

+ 0.005 in 

-1.02 

0.088 


0.0094 

-0.0012 

- 

+45 

-0.0029 


0.92 

in 

-0.0075 


(E7.5) 


(E7.6) 


(E7.7) 


(E7.8) 



Transforming the strains in the + 45° ply into the principal material directions. 


Top of ply No. 2: 



Top 

£ \ 


e 2 

= 

7n 


2 . 

+45 


cos 2 45° 
sin 2 45° 

- sin 45° cos 45° 


e^ 5 = 0.00029 ; 

Bottom of ply No. 2: 


sin 2 45° 
cos 2 45° 

sin 45° cos 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

|cos 2 45° - sin 2 45° j 


0.00425 

-0.00077 

-0.0029 

2 


e + 2 45 = 0.00319 ; = -0.00502 


(E7.9) 


* 

£ 1 

Bot 

£ 2 

= 

Yn 


2 

+45 


cos 2 45° 
sin 2 45° 

- sin 45° cos 45° 


gj 45 = 0.00035 ; 


sin 2 45° 
cos 2 45° 

sin 45° cos 45° 


2 sin 45° cos 45° 
-2 sin 45° cos 45° 

|cos 2 45° - sin 2 45°) 


0.0094 

- 0.0012 

-0.0075 

2 


e^ 5 = 0.00785 ; 7 ^ = -0.0106 


(E7.10) 


Putting these values of strain into tabular format: 


Strains 


Plane 

z (in) 

e x 

£y 

Yxy 

£\ 

e 2 

712 

Ply 1 (Top) 

+0.005 

-0.00085 

-0.00033 

0.0017 

-0.00085 

-0.00033 

0.0017 

Ply l(Bot.) 

0 

0.00425 

-0.00077 

-0.0029 

0.00425 

-0.00077 

-0.0029 

Ply 2 (Top) 

0 

0.00425 

-0.00077 

-0.0029 

0.00029 

0.00319 

-0.00502 

1 Ply 2 (Bot.) 

-0.005 

0.0094 

- 0.0012 

-0.0075 

0.00035 

0.00785 

-0.0106 


A diagram of the through-the-thickness strain in the x-direction follows: 




— — -s 





-1.0 -.8 -.6 --4 -.2 0 -2 .4 .6 .8 1.0 

£x (percent) 


(in.) 

.005 


0 


-.005 
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Note that although only a tensile stress is applied to the plate, a small amount of compressive strain 
is produced in the outer fibers of part of the 0° ply. This is due to the By terms which relate midplane 
strains to plate curvatures (i.e., due to the nonsymmetry, the plate is bending even though no 
bending moment is applied). As can be seen from the diagram, the plate is curving into the 0° fiber 
direction. 


Example 8: 

Consider an applied moment of 2.5 inch-lb/inch in the 0° fiber direction. The constitutive 
equation is: 


r 

0 


0 


0 


K x 


K y 




J 



^ 11 ^ 12^16 0 0 0 
^ 12 ^ 22^26 0 0 0 


^16^26^66 


0 0 0 


0 0 0 D n D i2 D \6 

0 0 0 ^ 12 ^ 22^26 


0 0 0 


D 16 D 26 D 66 


0 

0 

0 

M x 

0 

0 


From this the midplane strains and plate curvatures can be found. 


e * = -0.00407- 2.5 — 

* 11 * lb in 


- 0.0102 


( 


e =B, 0 M X = 0.000368 — 
y 12 lb 


2.5 


in-lb 


in 


= 0.00092 


r =B.,M X = 0.00375 

•xy 16 x 


lb 


2.5 — 1 = 0.00938 
in 


K X = D U M= 2.05— - 
11 in-lb 


K V = D,,M X = -0.366 — 
y 12 in-lb 


f ■ 

2.5 — — 




2.5 


V 


m 

in-lb 

in 


5.1 


in 


= -0.92 r 


in 


K =D. 6 M x = -\31 - — — 
^ 16 in-lb 


2.5 l _ _3.43 — 
in I m 


(E8.1) 


(E8.2) 
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Using equation (30), the strains through each of the two plies can be calculated. 

For the top of the 0° ply: 


£ . 

Top 

-0.0102 

0.00092 

+ 0.005 in 

1 

O 

KO t—k 
to 


0.0153 

-0.00368 


0 

0.00938 


-03.43 

in 

-0.0078 


(E8.3) 


For the bottom of the 0° ply (which will be the same as the top of the +45° ply since the laminate 
axes are being used, not the principal material axes for a specific ply): 


£ 

BoL 

-0.0102 

0.00092 

+ 0 in 

5.1 

-0.92 


-0.0102 

0.00092 


£ 

& - 

0 

0.00938 


-03.43 

in 

0.00938 


* 


(E8.4) 


For the bottom of the +45° ply: 


£ x 

£,, 

BoL 

-0.0102 

0.00092 

- 0.005 in 

5.1 

-0.92 


-0.0357 

0.00552 

y 

Yxy 

+45 

0.00938 


-03.43 

in 

0.0265 


(E8.5) 


Transforming the strains in the +45° ply into principal material directions: 


e i 

Top 

cos 2 45° 

sin 2 45° 

2 sin 45° cos 45° 

-0.0102 

£ 2 

= 

sin 2 45° 

cos 2 45° 

-2 sin 45° cos 45° 

0.00092 

Tl2 

2 

+45 

- sin 45° cos 45° 

sin 45° cos 45° 

cos 2 45° - sin 2 45°j 

0.00938 

2 



£^ 5 = 0.000 ; 

4-45 

£2 =-0-00933 

y£ = 0.0111 


£ 1 

Bot. 

cos 2 45° 

sin 2 45° 

2 sin 45° cos 45° 

-0.0357 

E 2 

= 

sin 2 45° 

cos 2 45° 

-2 sin 45° cos 45° 

0.00052 

CM 

+45 

- sin 45° cos 45° 

sin 45° cos 45° 

^cos 2 45°-sin 2 45°) 

0.0265 

2 


e + 45 = 0.00184 ; £^ 45 = —0.0283 ; y^ = 0.0412 


, (E8.6) 


. (E8.7) 
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Putting the above values of strain into tabular form: 


Strains 


Plane 

z (in.) 


£y 

Yxy 

£ i 

£ 2 

7\2 

Ply 1 (Top) 

+0.005 

1 

-0.00368 

-0.0078 

-0.00085 


0.0017 

Ply l(Bot.) 

0 

mm 

0.00092 

0.00938 

0.00425 


- 0.0029 

Ply 2 (Top) 

0 

-0.0102 

0.00092 

0.00938 

0 


0.0111 

Ply 2 (Bot.) 

-0.005 

-0.0357 

0.00552 

0.0265 

-0.00184 


0.0412 


A diagram of the through-the-thickness strain in the x-direction follows: 


Z (in.) 



.005 


0 


-.005 


Note that at the geometric midplane, the strain is not zero but -1.02 percent, even though 
only a bending moment has been applied. This is due to the bending-extension coupling that is 
present in this unsymmetric laminate. The plane of zero strain in the x-direction- can be determined 
by setting £ x equal to zero in equation(30) and solving for z: 


0 


-0.0102 


5.1 



0.00092 

0.00938 

+ z 

-0.92 

-3.43 


=> 0 = -0.0102 + |5.1 r^j z =>z= 0.002 in 


(E8.8) 


Thus, at z = 0.002 in, the x-direction strain is zero for this example. A diagram of the through-the- 
thickness strain along the fibers in each ply follows: 
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Z (in.) 



£l (percent) 


C. Determination of the Neutral Plane 

The plane of zero strain for any direction in an unsymmetric laminate can be calculated from 
equations (52) and (30). For pure bending in the x-direction: 


0 




D \\M X 

0 

= 

B n M x 

+ Z 

D n M x 

0 


/ 


/ 



B 16 ^x 


D 16 m x 


(96) 


B 


ll 


D 


ll 

r 

B 


12 


D 


12 


B 


16 


D 


16 


for the e = 0 plane , 


for the £y = 0 plane , 


for the y xy = 0 plane . 


(97) 
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For example No. 8, these values of z are; 


0.00407 


For thee =0 plane; z = 

* 2.05 


lb 


in-lb 

-0.000353 


For the e =0 plane ; z = 

y -0.366 


lb 


in-lb 
-0.00368 


For the r = 0 plane ; z = 

rxy v -L37 T 


lb 


= 0.002 in 


= 0.00096 in 


= 0.0027 in 


in-lb 


The general form of equation (96) is: 

For the £ x = 0 plane; 

0 = . ( 99 ) 


For the £y = 0 plane; 

' / f r / / 

0 = B^2M x +B22Mfy+B2^M X y+z(D^2^ x '^D22^y + D^^M X y)— 


in 


( 100 ) 


For the y xy = 0 plane; 

0 -B l6 M x +B2 6 M y +B 66 M xy +z(D l6 M x +D2 6 My+D 66 M xy )r^ . ( 101 ) 

It can be seen that if only one moment resultant is acting on the laminate then the planes of 
zero str s (the neutral planes) will not depend on the value. However, if more than one moment 
resultant is being applied to the plate, then the magnitude of these moment resultants will determine 
where the neutral planes are located. Regardless, equations (99), (100), and (101) can be used to 
find the planes of zero strain for an unsymmetric laminate subjected to bending moments. 


Example 9: 

Consider a square laminated plate of AS4/3501-6 with a stacking sequence of [+45, 0, 90]j\ 
Find the stresses and strains at room temperature (70 °F) in each ply in both the laminate and 
material directions due to thermal processing effects. The material was processed at 350 °F. 
Assume no moisture effects. 



The lamina properties are: 


E ! = 20,010,000 lb/in 2 , 
E 2 = 1,301,000 lb/in 2 , 
G n = 1,001,000 lb/in 2 , 
v 12 = 0.3 , 
v 21 = 0.02 , 
a x = -0.072x1 0 -6 — , 

a 2 = 32.4xl0 -6 — , 
Ply thickness = 0.005 in . 


Tom equation (10), the reduced stiffnesses can be calculated: 


Q 


it 


20,010,000 lb/m = 20439/785 lb/in 2 , 


!— Vi 2 v 2 i l-(0.3)(0.02) 


2 


Q 


22 ' 


E 2 1,301,000 lb/in" ., .. 2 

= 1,308,853 lb/m , 


1 “ v i2 v/ 2i 1 -(03)(0.02) 


v 10 £ 2 (0.3) (1,301,000 lb/in 2 ) 9 

Q 12 = — = - — - = 392,656 \b!m , 


l ~V\2 V 2\ 


l-(0.3)(0.02) 


2 66= 1 ,0° 1,000 lb/in . 


Tom equation (22); for the 0° plies: 


Q° n = 20,130,785 lb/in 2 , 
Q° n = 392,656 lb/in 2 , 
Q° 22 = 1,308,853 lb/in 2 , 

S°6 = o. 

526 = 0 , 

Q 0 66 = 1,001,000 lb/in 2 . 


(E9.1) 


(E9.2) 
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For the +45° ply: 


0 = (20,130,785 lb/in 2 )(0.707) 4 +2 


(392,656 lb/in 2 +2( 1,00 1,000 lb/in 2 ) (0.707) 2 (0.707)^ 


+(1,308,853 lb/in 2 )(0.707) 4 = 6,557,237 lb/in 2 , 


— +45 
012 “ 


r 2 


r 2 


(20,130,785 lb/in"+ 1,308,853 lb/in^-4( 1 ,00 1 ,000 lb/in 2 ) 1(0.707)^(0.707) 


/in 2 (( 


+392,656 lb/in" (0 4 +(0.707) 4 


)= 4,555,238 lb/in 2 , 


0|f = (20,130,785 lb/in 2 )(0.707) 4 +2^(392,656 lb/in 2 +2( 1,00 1,000 lb/in 2 ) j(0.707) 2 (0.707) 2 
+(1,308,853 lb/in 2 )(0.707) 4 = 6,557,237 lb/in 2 , 


— +45 


2 


0 l6 " = 20,130,785 lb/in"-392,656 lb/in 2 -2( 1,00 1,000 lb/in^) 


(0.707) (0.707) 


J 


+1 392,656 lb/in"- 1,308,853 lb/in"+2( 1,00 1,000 lb/in^) J(0.707)(0.707) J = 4,705,483 lb/in 2 , 


— +45 
&26 ~ 


20,130,785 lb/in 2 -392,656 lb/in 2 -2(l, 001,000 lb/in 2 ) 


(0.707) (0.707) 


J 


n.2 


r 2 


392,656 lb/in"- 1,308,853 lb/in"+2( 1,00 1,000 lb/in 2 ) |(0.707)(0.707) J = 4,705,483 Ib/in 


2 


Q* 45 

^66 


20,130,785 lb/in 2 + 1,308,853 lb/in 2 -2(392,656 lb/in 2 + 1,001,000 lb/in 2 ) 


(0.707) 2 (0.707) 2 


i-2 


For the 90° ply: 


+ 1,001,000 lb/in (0.707) +(0.707) = 5,163,582 lb/in 


gf = 1,308,853 lb/in 2 , 
0i2 = 392,656 lb/in 2 , 
022 = 20,130,785 lb/in 2 , 

qZ = 0 . 


• 16 ' 

226 = 0 . 


066= 1,001,000 lb/in 2 . 


(E9.3) 


n c 



From equation (53): 


An - 


20,130,785+6,557,237+1,308,853 


^ 12 “ 


A 22 - 


lb/in 2 (o.005 in) = 139,984 — , 
v ' in 

392,656+4, 555, 238+392, 656j lb/in 2 (o.005 in) = 26,703 — , 

1,308,853+6,557,237 + 20,130,785 lb/in 2 (o.005 in) = 139,984 — , 

v ' in 

A i 6 = [4,705,483] lb/in 2 (o.005 in) = 23,527 — , 

A 26 = [4,705,483] lb/in 2 (o.005 in) = 23,527 — , 


in 


^ 66 ~ 


1,001,000+5,163,582+1,001,000 


lb/in 2 (0.005 in) = 35,828 — 
v ' in 


From equation (54); the Bi/s will be: 


B n= [6,557,237(0.005)— 1,308,853(0.005)] ^-(0.005 in) = 131.21b , 
B i 2 = [4,555, 238(0.005)-392,656(0.005)] jjj- (0.005 in) =104 lb , 

B 22 = [6,557, 237(0.005)-20, 130,785(0.005)] ^-(0.005 in) =-339.4 lb , 

B j 6 = [4,705,483] lb/in 2 (0.005 in) (0.005 in) = 117.61b , 

B 26 = [4,705,483] lb/in 2 (0.005 in) (0.005 in) = 117.6 lb , 

B 66 = [5,163,582(0.005)-1,001,000(0.005)] ^-(0.005 in) = 1041b . 


From equation (56); the £>, -/s will be: 


(E9.5) 


(E9.6) 
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D 


ll 



( 3 1 


r 3 ) 

6,557,237 

— — + (0.005) (0.005) 2 

+ 20,130,785 

(0.005) 3 +q 


l 12 J 


, 12 j 


( . i 


\ 


+ 1,308,853 


V 


( °'^ 5) - + (0.005) (0.005) 2 


J 


1.275 in- 


£> 12 “ 


4,555,238 


V 


( °'^ 5) + (0.005) (0.005) 2 | + 392,656^ ( °'^ 5) + 0 
(0.005) + (0 .005) (-0 .005) 2 


+ 392,656! 


V 


12 


D 22~ 


( 


6,557,237i 




(0 - 005) + (0.005) (0.005) 2 


V 


12 


+ 1,308,853 


J 


f (0.005) 3 


= 0.674 in-lb , 


J 




12 


+ 0 


! 


+ 20,130,785 


V 


(0 °° 5) + (0.005) (-0.005) 2 


) 


= 3.628 in-lb , 


£>16 = 


f 


4,705,483 


3 


V 


- °' - ^ 5) + (0.005) (0.005) 2 


J 


= 0.637 in-lb , 


D 26 = 


4,705,483 


(°- 005) + (0.005) (0.005) 2 


V 


12 


J 


= 0.637 in-lb , 


£> 66 “ 


5,163,5821 


(0 - 005) + (0.005) (0.005) 2 


V 


12 


1 , 001,000 


f 


(0.005)- 


1 , 001,4 (Q - QQ5) +0 

V 12 
3 

2 


V 


12 


+ (0.005) (-0.005) 


J 


■ 0.845 in- 



The extensional stiffness [A], coupling stiffness [B], and bending stiffness [£>] matrices are: 

139,984 26,703 23,527 
26,703 139,984 23,527 
23,527 23,527 35,828 


131.2 104 117.6 

[b] = 104 -339.4 117.6 lb , 

117.6 117.6 104 


1.275 0.674 0.637 
0.674 3.628 0.637 in-lb . 

0.637 0.637 0.845 

The fully inverted form of the constitutive equations is needed. First, the partially inverted form is 

given by equation (51): 

8.1 -0.73 -4.8 
= [a] _1 = -0.73 8.1 -4.8 
-4.8 -4.8 34.3 

-4.22 -5.26 -3.67 
-1.82 33.9 -3.67 
-29.05 -51.6 -24.38 


4.22 1.82 29.05 

[C'U*M 1 = 5.26 -33.9 51.6 xlO^in , 
3.67 3.67 24.38 

0.859 0.350 0.264 

[4M-WW1 = 0.350 1.816 0.437 in-lb . 

0.264 0.437 0.505 
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The fully inverted form can now be determined from equation (52a): 




8.5 

-0.91 

-2.175 

-0.91 

18.1 

-10.3 

-2.175 

-10.3 

55.2 


1 A" 6 111 

x 10 — 
lb 




* . 

D = 


-299 -120 -467 
-427 2,614 -2,765 
-1,932-1,959-2,124 


x 10 


-6 


lb 


[c']~[aTV]= 


-299 -427 -1,932 
-120 2,614 -1,959 
-467 -2,764-2,124 


x 10 


-6 


lb 


(E9.8) 



1.41 -0.119-0.633 
-0.119 0.706 -0.548 
-0.633-0.548 2.79 


in-lb 


The thermal expansion coefficients for each ply can be calculated from equation (90): 
For the +45° ply: 

a i =16.16xl(r 6 -: a y = 16 . 16 xl 0 -6 — ; ^ = -16.24xlCf 6 F . 

For the 0° ply: 

a = -0.072xl0 -6 — ; cc = 32.4xl0~ 6 — ; a =0— . 

x p° y p° xy po 

For the 90° ply: 

^ = 32.4xl0- 6 - ; ^ = -0.072x10-*- ; % = 0- . 

The thermal stresses in each ply can be found using equation (91). 
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For the +45° ply: 


e=a AT = 16.16x10 6 — (-280 °f) = -0.004525 , 
e y = cc y AT = 16. 16x 10" 6 — (-280 °f) = -0.004525 , 
y = 2a Ar = -32.48xl0” 6 — (-280 °f) = 0.009094 . 

'xy xy p° v 1 

For the 0° ply: 

/ = a x AT = -0.072xl0~ 6 — (-280 °f) = 0.00002016 , 
/ = a y AT= 32.4xl0 -6 — (-280 °f) = -0.009072 , 

Y =2a AT=0 — (-280°f) = 0 . 

'xy xy po V 1 

For the 90° ply; 

e x = a x AT= 32.4xl0“ 6 — (-280 °c) = -0.009072 , 
/= a AT = -0.072xl0" 6 — (-280 °c) = 0.00002016 , 

y y v— 

i = 2 «^7-=0-(-280» c ) = 0 . 

The thermal stress resultants can be found using equation (92): 


(E9.9) 


(E9.10) 


(E9.ll) 
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T 

"x 



T 

N y 

II 


3 ^ 




6.557 4.555 4.705 
4.555 6.557 4.705 
4.705 4.705 5.164 

+ 

20.131 0.393 0 

0.393 1.309 0 

0 0 1.001 

1.309 0.393 0 

0.393 20.131 0 

0 0 1.001 


-0.004525 

-0.004525 

0.009094 


0.00002016 

-0.009072 

0 

-0.009072 

0.00002016 

0 


+ | xlO 6 (0.005 in) lb/in 2 


- 112.6 

- 112.6 

21.91 


lb 

in 


(E9.12) 


The thermal moment resultant is given by equation (93): 


V 



T 

M y 

= 


T 

KJ 


r 


6.557 4.555 4.705 
4.555 6.557 4.705 
4.705 4.705 5.164 


-0.004525 

-0.004525 

0.009094 


0.1093 

-0.1084 

0.1095 


lb 


(o.005 in) 


20.131 0.393 0 

0.393 1.309 0 

0 0 1.001 

1.309 0.393 0 

0.393 20.131 0 

0 0 1.001 


0.00002016 

-0.009072 

0 


(o in) 


-0.009072 

0.00002016 

0 


(-0.005 in) 


xlO 6 (0.005 in) lb/in 2 
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The mechanical stress and moment resultants have only one component due to the applied load; M x 
= -16.7 lb. Thus, the complete set of stress and moment resultants are: 


N. 


Tot 


N 


X 

Tot 


N. 

m: 


y 

Tot 


M. 


xy 

Tot 

x 

Tot 


M 


y 

Tot 


xy 


- 112.6 

- 112.6 

21.91 


lb 

in 

lb 

in 

lb 

in 


0.10931b 
-0.1084 lb 
0.10951b 


(E9.14) 


The midplane strains and curvatures can now be calculated from the constitutive equations (52). 

£ ° x = [(8.5X-1 12.6) - (0.91X-1 12.6) - (2.175)(21.91)) xl0~ 6 

+ [(-299)(0. 1093) -(120X-0. 1084) -(467)(0. 1095)] xlO -6 = -0.0009731 

e = [(— 0.9 1)(— 1 12.6) + (18.1X-112.6)- (10.3X21.91)] xlO -6 

+ [(-427)(0.1093) + (2,614)(-0.1084) - (2,765)(0.1095)] xlO -6 = -0.002794 

Y .= [(-2. 175)(-1 12.6) - (10.3)(— 1 12.6) + (55.2)(21.91)] xlO -6 

+ [(— 1,932)(0.1093) - (1,959)(-0.1084) - (2,124)(0.1095)] xlO -6 = 0.002383 

(E9.15) 

K x = [(— 299)(— 1 12.6) - (427)(-l 12.6) - (1,932)(21.91)] xlO -6 

+ [( 1 .4 1 )(0. 1093) - (0. 1 1 9)(-0. 1 084) - (0.633)(0. 1095)] = 0.1371 — ’ 

m 

K y = [(— 120)(— 1 12.6) + (2,614)(-1 12.6) - (1,959)(21.91)] xlO -6 

+ [(-0. 1 19)(0. 1093) + (0.706)(-0. 1084) - (0.548)(0. 1095) = -D.4733 — ’ 
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= [( — 467)(— 1 12.6) - (2,764)(-l 12.6) - (2,124)(21.91)] xlO 6 

+ [(-0.633)(0. 1 093) - (0.548)(-0. 1 084) + (2.79X0.1095)] = 0.6130 - 

in 

Now equation (30) can be used to find the strains within the laminate, and once these are 
found, equation (23) can be used to find the stresses. 

At the top of the +45° ply, z = 0.0075 in: 


e * 

£ v 

Top 

-0.0009731 

-0.002794 

+ 0.0075 in 

0.1371 

-0.4733 


0.000055 

-0.006344 

y 

445 

0.002383 


0.6130 

in 

0.006981 


At the bottom of the +45° ply (the top of the 0° ply, z = 0.0025 in: 


£ v 

Bot 

-0.0009731 

-0.002794 

+ 0.0025 in 

0.1371 

-0.4733 


-0.000630 

-0.003977 

y 

445 

0.002383 


0.6130 

in 

0.00391 


At the bottom of the 0° ply (the top of the 90° ply), z = -0.0025 in: 



Bot 

-0.0009731 


0.1371 


-0.001316 

£ y 

0 

-0.002794 

0.002383 

- 0.0025 in 

-0.4733 

0.6130 

in 

-0.001611 

0.000851 


At the bottom of the 90° ply, z = -0.0075 in: 


£ v 

Bot 

-0.0009731 

-0.002794 

- 0.0075 in 

0.1371 

-0.4733 


-0.002001 

0.000756 

y 

90 

0.002383 


0.6130 

in 

-0.002215 


(E9. 16) 


(E9.17) 


(E9.18) 


(E9.19) 
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Transforming the +45° ply strains into principal material strains: 


£ 1 

Top 

0.5 0.5 1 


0.000055 


0.000346 

£ 2 

n 2 

2 

+45 

0.5 0.5 -1 
-0.5 0.5 0 


-0.006344 

0.006981 

2 


-0.006635 

-0.003200 


(E9.20) 


£ i 


0.000346 

e 2 

= 

-0.006635 

ru 


-0.006400 




£ i 

But. 

0.5 0.5 1 


-0.00063 


-0.000346 

£ 2 

rn 

2 

+45 

0.5 0.5 -1 
-0.5 05 0 


-0.003977 

0.003916 

2 


-0.004262 

-0.001674 


(E9.21) 



£ i 


-0.000346 

=> 

£ 2 

= 

-0.004262 


Yll 


-0.003348 


The stresses in the laminate material directions are; 

Note: Recall that the thermal strains must be subtracted out of the total strain. 


F 1 

°x 

Top 





1 



6.557 4.555 4.705 


0.000055-(-0.004525) 


11.80 


G y 

T xy 

+45 

4.555 6.557 4.705 
4.705 4.705 5.164 

xlO 6 lb/in 2 

-0. 006344-(-0. 0045 25 ) 
0.006981-0.009094 


p g 
7 

ksi , (E9.22) 

°x 

Bot 

6.557 4.555 4.705 


-0.000630-(-43. 004525) 


3,67 


T 


4.555 6.557 4.705 
4.705 4.705 5.164 

xlO 6 lb/in 2 

-0.003977-(-0.004525) 

0.003916-0.009094 


-3.03 

-5.83 

ksi , (E9.23) 

V 

xy 

+45 

- 


- 





■ 

Top 

- 


- 


■ 


■ 



20.131 

0.393 

0 


-0.000630-0.00002016 


-11.09 

°y 

= 

0.393 

1.309 

0 

xlO 6 lb/in 2 

-0.003977-(-0. 00907 2) 

= 

6.41 

T 


0 

0 

1.001 


0.003916-0 


3.92 

Z xy 

0 

- 


■ 


• 


l J 


■ 

Bot 

- 




r -I 


r 

a x 


20.131 

0.393 

0 


-0.001316-0.00002016 


-24.0 


= 

0.393 

1.309 

0 

xlO 6 lb/in 2 

-0.00161 1-(-0.009072) 

= 

9.24 

T 


0 

0 

1.001 


0.000851-0 


0.85 

Z xy 

0 

- 


• 


• 


L 


* 

Top 

- 


■ 


■ 


r i 

ff x 


1.309 

0.393 

0 


-0. 00 1 3 1 6-(-0. 009 072) 


9.51 

s 

= 

0.393 

20.131 

0 

xlO 6 lb/in 2 

-0.001611-0.00002016 

= 

-29.8 

T 


0 

0 

1.001 


0.000851-0 


0.85 

T xy 

90 

- 


■ 


- 


L J 


■ 

Top 

- 


- 


■ 


r 

°x 


1.309 

0.393 

0 


-0.002001-(-0.009072) 


9.55 

°y 

= 

0.393 

20.131 

0 

xlO 6 lb/in 2 

-0.000756-0.00002016 

= 

17.6 

T 


0 

0 

1.001 


0.002215-0 


-2.22 

xy 

90 

■ 


• 


• 


L J 


Transforming these stresses into the +45° ply’s principal material directions; 


a i 

Top 


= 

T 12 

+45 


0.5 

0.5 

-0.5 


0.5 

0.5 

0.5 


11.80 


7.48 


-1.01 

ksi = 

3.32 

ksi , 

2.08 


-6.41 




Bot 


= 

T n 

+45 



0.5 

0.5 

-0.5 


0.5 

0.5 

0.5 


3.67 


-5.51 


-3.03 

ksi = 

6.15 

ksi 

-5.83 


-3.35 



(E9.24) 


(E9.25) 


(E9.26) 


(E9.27) 


(E9.28) 


(E9.29) 


Putting the strain and stress results into tabular form: 


Ply Strains 


Location 

e x 

£y 

Jxy 

€l 

£ 2 

Yl2 

Top +45 
Bot. +45 

0.000055 

-0.006344 

0.006981 

0.000346 

-0.006635 

-0.006400 

-0.000630 

-0.003977 

0.003916 

-0.000346 

-0.004262 

-0.003348 

TopO 
Bot. 0 

-0.000630 

-0.003977 

0.003916 

-0.000630 

-0.003977 

0.003916 

-0.001316 

-0.001611 

0.000851 

-0.001316 

-0.001611 

0.000851 

Top 90 
Bot. 90 

-0.001316 

-0.001611 

0.000851 

-0.001611 

-0.001316 

-0.000851 

-0.002001 

0.000756 

-0.002215 

0.000756 

-0.002001 

0.002215 

Ply Stresses in ksi 

Location 

Ox 

C7 y 

Try 

CXi 

02 

Tl2 



-1.01 

2.08 

7.48 

3.32 

-6.41 

Bot. +45 

3.67 

-3.03 

-5.83 

-5.51 

6.15 

-3.35 



6.41 

3.92 

-11.09 

6.41 

3.92 

■ 

■ ■ 

9.24 

0.85 

-24.0 

9.24 

0.85 

E 


9.51 

-29.8 

0.85 

-29.8 

9.51 

-0.85 

9.55 

17.6 

-2.22 

17.6 

9.55 

2.22 

Bot. 90 


The deformed shape of the plate can be predicted from the midplane strain and plate curvature 

results: 

e° x = -0.0009731 , 

£° = -0.002794 , 
r° = 0.002383 , 

' xy 

K x = 0.1371 , (E9.30) 

£ =-0.4733- , 
y in 

£ =0.6130- . 
y in 


From the definition of plate curvatures, equation (29), it is seen that these can be used to find the 

out-of-plane displacements w. 

Along the edge parallel to the 0° fibers (the x-direction) the plate will bend with a change in 

w of: 

w(x) = J 2 J “-£* = -0.1371^ (* 2 -*i] + vv o(-* 2 _x i) + w o • (E9.31) 
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Let wo = 0 and xi = 0. This edge will vary parabolically with length as shown in the following: 

x (in.) 



Along the edge perpendicular to the 0° ply fibers (the y-direction): 


w(y) =JT 2 2 -^ = -°- 4733 ^(y 2 ->’i ) 2 + w o(> ? 2-) ; i) + vt; o (E9.32) 

This edge will also vary parabolically, but with a different bending direction with a more pronounced 
curvature. 


w (y) in. 
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There will also be a twist to the plate: 


w(x,y) = j^ 2 \ y 2 = -0.613 - (y 2 -y j)^-* i) + vv 0 (* 2 -* i) + w o • (E9.33) 

Again letting wo = 0 and x\ and yi = 0: 

w = 3.7 in. 



1 2 3 


Superimposing the three curvatures will demonstrate what the plate will look like after it cures to 

room temperature: 



IX. REFERENCE OF IMPORTANT EQUATIONS 


Transformation matrix: 


M- 


2 Z1 

cos 9 
sm 9 


sin 

2 ^ 
cos 9 


sin 6? cos 9 sin 0cos 9 


2 sin 9 cos 9 

-2 sin 9 cos 9 

2 „ . 2 . 

cos 9 - sm 9 


< 7 . 




H"= 


2 n 

cos 9 
2 

sm 9 


■ 2 Zi 

sin 9 
2 

cos 9 


- 2 sin 0cos 8 
2 sin 9 cos 9 


2 - 2 

sin0cos0 -sin0cos0 Icos 0-sin 6 


<*l 


°x 

°2 

= m 

c y 

►— * 
to 


T xy 


<*x 

= [T\~ l 

°i 

°y 

®2 

£ 


hi 


h 



£ 1 

= m 


S 12 





= [71" 1 

e i 
& 2 

$ 


£ 12 
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Stress-strain relationship for a lamina, along the principal material directions: 


°l 


QuQn o 

e l 

a 2 

= 

2 12 222 0 

£ 2 

T 12 


o o 266 

n 2 


Qi 1 = 


l ~ y n v 2 1 


2 


22 ' 


l -V\2 y l\ 


Qn~ 


V 12 E 2 


V 21 E l 


l-v n v 2l l-v 12 v 21 


266-^12 


Stress-strain relationship for a lamina: 


Qu = Qn m4+2 (Qn +2 Q66) m2n2+ Q22 n4 
Ql2=[Qn + Q22- 4 Q66] m2n2+ Qu( m4+n4 ) 
Q22 = Qll n + 2^Q 12 + 2Q 6e] m n + Q22 m 
Ql6~ {Qn~Q 12~^Q 66 n+ {Q 12~2 22+22 66 j mn 
Q26 = [Q 11~2 12 — ^2 66 m+ ^Q 12“2 22+^2 6e) nm 
Q(,6~ [Qll + Q22~2Ql2~2Q6<^ m n + 26e( m +n ) 


m = COS 6 

* 

n = sin 6 


°x 


2 11 2 12 2l6 



a y 

II 

2 12 2 22 2 26 



T xy 


2 16 226 266 
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Laminate constitutive equations: 




1 



An A 12 A 16 | #11 # 12 



A 12 A 22 A 26 1 #12 #22 



^16 ^26^66 J # 16 B 26 



#11 # 12 # 16 | #*11#*12 



#12 #22 #26 1 #*12 #*22 



# 16 #26 #66 j #*16 #*26 


B 

B 

B 


D 

D 

D 


16 

26 

66 


16 

26 

66 


K , 


K 




Extensional stiffness matrix: 


Coupling stiffness matrix: 


Bending stiffness matrix: 


Partially inverted forms: 





kl-w 1 . 

k]=-W'M . 


c 





OK 



Fully inverted forms: 


/ 


4» 



A 

-l 


A 

= 

*T* 

A 

- 

•T* 

5 

T 

D 


C 


/ 


* 

* 

B 

= 

5 

D 


f 


* 

-1 

* 

C 

zz - 

D 


c 


/ 


* 

D 

= 

D 


Fully inverted constitutive equations: 






An Aj 2 Aj 6 
/ / / 

A 12 A 22 ^26 

/ / r 

A 16 A 26 A 66 


5 11 5 12 B 16 
/ t t 

B \2 B 22 B 26 
/ / / 

B 16 B 26 B 66 


Ny 

Hxy 

M 
M y 
M 


K > 

K 


xy 


C 11 C 12 C 16 
C \2 C 22 C 26 


D U D 12 D 16 
' / • 

D \2 D 22 D 26 

Ci 6 C 26 C 66 1 d\ 6 d' 26 d' 66 


xy 


Strains within a laminate: 


i 


o 







K x 

£ 

=: 

o 

£ 

+ z 

K v 

y 


y 


y 

r*y 


o 

7 

'xy 


a 


Engineering constants for symmetric laminates: 


Ah A 12 


h 


A 26 A 16~ A n A 




12^66 


A 22 A 66~ A 26 


+ ' 


16 


AofiAioAf 


16 ^26^ 12^66' 
+ - 


- A 26 A 16 


1 66 


A 22 A 66 


~ A 26 A 66 


17 ^ 12 


! 


A i^A 2 g— A 1 7 A 


12^66 


^ 1 1 ^ 66 “^ 


16 


+ ' 


‘22 


+ ' 


1 26 


l 26 


+ 


■66 


A 16 A 12 A 66~~ A 16 A 26 
2 2 ‘ 

A 1 1^66 ~ A 16 ^66 


a 


92 



A 66 A 26 2A n A 16-^26^22"^ 12 A 26 ~A 16 A 22 
h W A 22 h (a 1 ! A 22 2 -A j 2 2 A 22 ! 


Thermal stains: 


*y 


^ 12 “ 


A 16^26 
A 66 


f 

A 22 

V 



A 66 J 


yx 


A 16 A 26 

^ A 66 


\ 

12 

) 



\ 


e x = ot x AT , 


cc y AT 


/ = — = a AT 

^2 ** 


Thermal stress and moment resultants: 


r 

n 

2n 2i22i6 

r 

Ny 

-2 

2 12 022 026 

n t 

/t=i 

2 16 226 266 

xy 


• 



n 

2 11 2i 2 2i6 


< 

-2 

2 12 222 226 


T 

M 

*:=i 

2l6 226 266 


xy 


■ 

k 


r. 


*y 




y 

r 


r 5 


*y 


k 
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Total stress and moment resultants: 


nT 


T 

N x 


n x 


f ju Tot 
M x 


T 

M x 


M x 

Tot 

N 

y 

= 

T 

N y 

+ 

N J 

and 

S/f Tot 

M 

y 

= 

T 

M y 

4 - 

M y 

N Jot 

xy 


T 

In 


N *y 


4 £ 


T 

rA 


M *y 


94 



REFERENCES 


1. Jones, R.M.: “Mechanics of Composite Materials.” New York: McGraw-Hill, 1975. 

2. Halpin, J.C.: “Primer on Composite Materials: Analysis.” Revised edition, Technomic 
Publishing Co., Inc., 1984. 



APPENDIX 


A Fortran program to find the determinant of a 6x6 matrix is given below: 

C THIS PROGRAM CALCULATES THE DETERMINANT 
C OF A 6 X 6 MATRIX WITH ELEMENTS "E" 

C 

REAL E(6,6), ES(6,6), M 
INTEGER TAG,I,J,N 

M=1 

N=1 

DO 101=1,6 

PRINT *, 'ENTER ROW #', I, 'OF THE MATRIX' 

READ *, E(1, 1 ),E(I,2),E(I,3),E(I,4),E(I,5),E(I,6) 

10 CONTINUE 
40 CONTINUE 
TAG=N 
DO 75 I=N,5 

IF (ABS(E(TAG,N)) .GT. E(I+1,N)) THEN 

TAG=TAG 
ELSE 
TAG =1+1 
ENDIF 

75 CONTINUE 
DO 80 J=N,6 
ES(N,J)=E(TAG,J) 

ES(TAG,J)=E(N,J) 

80 CONTINUE 
DO 90 J=N,6 
E(N,J)=ES(N,J0 
E(TAG,J)=ES(TAG,J) 

90 CONTINUE 
M=M*E(N,N) 

DO 91 J=N,6 

ES(N,J)=E(N,J)/E(N,N) 

91 CONTINUE 
DO 95 J=N,6 

E(N,J)=ES(N,J) 

95 CONTINUE 

DO 200 I=N+1,6 
DO 150 J=N,6 

ES(I,J)=(- 1 *E(I,N)*E(N,J))+E(I,J) 

150 CONTINUE 
200 CONTINUE 

DO 300 I=N+1,6 
DO 250 J=N,6 
E(I,J)=ES(I,J) 

250 CONTINUE 
300 CONTINUE 
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N=N+1 

IF (N .LT. 6) THEN 
GO TO 40 
ELSE 

CONTINUE 

ENDIF 

M=M*E(6,6) 

PRINT *, 'DETERMINANT =', M 

STOP 

END 


The determinant of a 5x5 matrix can be determined in a similar manner, substituting 5 for 6 in the 
appropriate sections of the above program. 
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